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Abstract 

The lectures are a practical introduction to perturbative calculations in QED 
and QCD. I discuss methods of calculation of one- and two-loop diagrams in di- 
mensional regularization, MS and on-shell renormalization schemes, decoupling of 
heavy-particle loops. 



Contents 



ll Introducti 



ion 



2 One-loop diagrams 



2.1 Massive vacuum diagram 



2.2 Massless propagator diagram 



3 QED at one loonl 



3.1 Lagrangian and Feynman rulesi 

3.2 Photon propagatod 

3.3 Ward identitvl 



3.4 Photon self-energv 



3.5 Photon field rcnormalizationi 



3.6 Electron propagator 



3.7 Vertex and charge renormalization 



3.8 Electron mass 



4 

4 
7 

10 

10 

11 

13 
15 
16 
18 
21 
26 



^ Given at Dubna International Advanced School on Theoretical Physics, 29 Jan. - 6 Feb. 2005 



1 



4 QCD at one loop 



4.1 Lagrangian and Fevnman rulesi 



4.2 Quark propagatoi 



4.3 Gluon propagatoi 



4.4 Ghost propagatoi 



4.5 Quark-gluon vertex 



4.6 Coupling constant renormalization 



4.7 Ghost-gluon vertex! 



5.2 Photon self-energvl 



Two-loop corrections in QED and qcdI 



5.1 Massless propagator diagram 



5.3 Photon field renormalizationi 



5.4 
5.5 



Charge renormalization 
Electron self-energvl . . 



5.6 Electron field renormalizat 

5.7 Two-loop corrections in QC 



iorJ 

53 



6.1 On-shell renormalization of photon field 



6.2 One-loop massive on-shell propagator diagram 



On-shell renormalization schemel 



6.3 On-shell renorrnalization of electron mass and fieldl 

6.4 On-shell chargel 



6.5 Magnetic moment 



6.6 Two- loop massive vacuum diagram 



6.7 On-shell renormalization of photon field and charge at two loopsi 

6.8 On-shell renormalization in OCDI 



Decoupling of heavv-particle loopsi 



7.1 Photon field 



7.2 Electron fieldl 



7.3 Chargel 



7.4 Electron mass 



7.5 Decoupling in OCD 



8 Conclusion: Effective field theoriei 

A Colour factord 

Referenced 



1 Introduction 



These lectures are an addendum to a standard textbook on quantum field theory, such 
as pp. We shall discuss practical ways to calculate one- and two- loop diagrams in QED 
and QCD and to obtain renormalized physical results. More details about calculation of 
Feynman integrals can be found in my lectures at another Dubna school 0, or in the 
textbook P^. Calculation of colour factors of diagrams in QCD (Appendix |XJ , and group 
theory in general, is discussed in an excellent book by P. Cvitanovic |lj available on the 
Web. In Sect. 13 about decoupling of heavy-particle loops, we follow [Sj, but in a much 
simpler case (QED at two loops). 

It is well-known that Feynman diagrams with loops can diverge. The first thing we 
have to do is to introduce some regularization. Regularized diagrams must converge, so 
that working with them has sense. When the regularization parameter tends to some 
value, the original theory must be restored. There exist many regularization methods. A 
good choice of regularization is crucial for any non-trivial computation: the regularization 
should preserve as much of the symmetry of the theory as possible, manipulating regular- 
ized Feynman integrals should be simple, etc. Cut-off regularization is difficult to define 
unambiguously; in gauge theories, it breaks gauge invariance; it is difficult to integrate 
by parts because of boundary terms. Pauli-Villars regularization also breaks gauge invari- 
ance. Lattice regularization preserves it, but breaks Lorentz invariance. The method used 
in practically all multiloop calculations is dimensional regularization. It preserves both 
Lorentz and gauge invariance, and rules for manipulating integrals are simple (no bound- 
ary terms, etc.). We shall use it throughout these lectures. We'll discuss the most widely 
used method of multiloop calculations, integration by parts [Hj, in Sect. 15. ll 

Another important concept is renormalization. Physical results should be expressed via 
physical parameters (measurable, or expressible via measurable quantities), not via param- 
eters in the bare Lagrangian. This re-expressing in physically necessary, independently of 
presence or absence of divergences. After it, we can switch off regularization; in a sensible 
theory, results must be finite. 

Renormalization procedure is not unique. There are many renormalization schemes; 
some of them depend on parameters. The most widely used scheme is MS. It depends 
on a single parameter — renormalization scale. Physical results should not depend on 
this artificial quantity. This requirement leads to renormalization group equations — a 
very powerful instrument in quantum field theory. We shall use this scheme in Sects. OHHl 
On-shell renormalization scheme (Sect. EI) is also widely used, especially in QED at low 
energies. In QCD, it is often useful to use on-shell renormalization for heavy quarks, and 
MS for all the rest. 
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2 One-loop diagrams 



2.1 Massive vacuum diagram 

During these lectures, we are going to live in dimensional space-time: 1 time and d — 1 
space dimensions. The dimensionality d must appear in all formulas as a symbol, it is not 
enough to obtain separate results for a few integer values of d. After calculating a physical 
result in terms of physical parameters, we take the limit c? — > 4. Therefore, d is often 
written as 4 — 26, where the physical limit is £ ^ 0. It is not possible to take this limit 
in intermediate formulas, because they often contain poles 1/e. Such divergences must 
disappear in physical (renormalized) results. 



The power of m is evident from the dimensional counting, and our aim is to find the 
dimensionless function V{n)] we can put m = 1 to simplify the calculation. Due to the 
position of the poles in the complex ko plane (determined by —iO), we may rotate the 
integration contour counterclockwise by 7r/2 without crossing the poles. After this Wick 
rotation, we integrate along the imaginary axis in the /cq plane in the positive direction: 
kf) = ikf). Here feo is the 0-th component of the vector k in d-dimensional euclidean 
space (euclidean vectors will be denoted by the bold font), and A;^ = —k^. Then our 
definition (j2.H) of V{n) becomes 



It is often useful to turn denominators into exponentials using the a-parametrization 



O 



Figure 1: One-loop massive vacuum diagram 



Let's consider the simplest diagram shown in Fig. ^ 




(2.1) 




(2.2) 




(2.3) 



For our integral, this gives 
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The d-dimensional integral of the exponent of the quadratic form is the product of d one- 
dimensional integrals: 



(2.5) 



This is the definition of (i-dimensional integration; note that the result contains c? as a 
symbol. Now it is easy to calculate 



V{n) 



T n) 



the result is 



V{n) 



T{-d/2 + n) 



For all integer n, the results are proportional to 

4 



Vi = Vil) 



{d-2){d-A) 



ra + e) 



(2.6) 



(2.7) 



(2.^ 



where the coefficients are rational functions of d. 

The denominator in ()2.2|) behaves as (fc^)" at k ^ oo. Therefore, the integral diverges 
if > 2n. At — s> 4 this means n < 2. This ultraviolet divergence shows itself as a 1/e 
pole in (j2I7|) for n = 1, 2. 

Let's calculate the full solid angle in d-dimensional space (we shall need it later). To 
this end, we'll calculate one and the same integral in Cartesian coordinates 



e '^''dk^ 



TC 



d/2 



and in spherical coordinates 
Therefore, the full solid angle is 



e-^'k'^-Mk 



^d 
2 



e-''\ef^-'de 

27r^/2 



ndT{d/2) 



(2.9) 



(2.10) 



(2.11) 
(2.12) 



T{d/2) 

For example, 

ni = 2, ^2 = 27r , ^3 = 47r , ^4 = 2??^ ... 

In 1-dimensional space, a sphere consists of 2 points, hence Qi = 2; the values for d = 2 
and 3 are also well known. 

What is the value of the massless vacuum diagram 

d'^k 



5 




Figure 2: One-loop massless vacuum diagram 

(Fig. 12))? Its dimensionality is d — 2n. But it contains no dimensionfull parameters from 
which such a value could be constructed. The only result we can write for this diagram is 

0. (2.13) 



This argument fails at n = d/2; surprises can be expected at this point. 




Figure 3: A diagram for the quark propagator 

This dimensions-counting argument is much more general. For example, the diagram 
in Fig. 01 contains a sub-diagram which is attached to the rest of the diagram at a single 
vertex, and which contains no scale. This subdiagram is given by an integral (maybe, a 
tensor one) with some d-dependent dimensionality. It has no dimensionfull parameters; 
the only value we can construct for such an integral is 0. 

Now we shall discuss a trick (invented by Feynman) to combine several denominators 
into a single one at the price of introducing some integrations. We'll not use it during 
these lectures; however, it is widely used, and it is important to know it. Let's multiply 
two a-parametrizations ()2.3j] : 



1 



a'l'a^' r(ni)r(n2) 



-aiai—a2a2 'ni — l n2 — l 



a'l^''afdaida2. (2.14) 



It is always possible to calculate one integral in such a product, namely, to integrate in a 
common scale rj of all a^. If we denote t] = ai + a^-, i.e., make the substitution a\ = r]x, 
a2 = ri{l — x), then we arrive at the Feynman parametrization 



T(ni + n2) x"i-^(l -x)"2-ic/x 



a^a^ T{n,)r{n2) Jo [a^x + a^il - x)]"^+"^ 



(2.15) 
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This is not the only possibihty. We can also take, say, t] = a2, i.e., make the substitution 
«! = rjx, a2 = Tj. This results in a variant of the Feynman parametrization, which can be 
more useful in some cases: 



Let's consider the general case of k denominators: 

1 _ 1 

a^^---a^'= ~ r(ni)r(n2) ■ ■ ■ T^) 



(2.16) 



(2.17) 



k ■ 



We can choose the common scale r] to be the sum of any subset of af, the numbering 
of the denominators is not fixed, and we can always re-number them in such a way that 
rj = ai + a2 ■ ■ ■ + <yi, where 1 < / < A;. We insert 

6{ai + a2 h a« - r])dr] 

into the integrand, and make the substitution ctj = r]Xi. Then the integral in t] can be 
easily calculated, and we obtain the general Feynman parametrization 

1 1(^1+^2 hnfc) 



aTar---al^ r{n,)r{n2) ■ ■ -Tin,) 

5{xi + X2 ■ ■ ■ + Xi — l)x"^~^X2^~^ ■ ■ ■ dxi dx2 ■ ■ ■ dXk 



X 



[aiXi + a2X2 h OfcXfc 



Let's stress once more that the 6 function here can contain the sum of any subset of the 
variables Xj, from a single variable up to all of them. 

2.2 Massless propagator diagram 

Now we shall consider the massless propagator diagram (Fig. 0}, 

/d^ = ^^'^'(~^'^'^''"^^^^^'^^^' ^i = -(^ + p)'' ^2 = -A;^ (2.19) 

(from now on, we'll not write —iO in denominators explicitly, but they are implied). We 
shall consider this diagram at p"^ < 0, i.e., below the threshold of production of a real pair, 
where the result is an analytic function. The power of —p'^ is evident from dimensionality. 
Our aim is to calculate the dimensionless function G{ni,n2)] we can put — = 1 to 
simplify things. This diagram is symmetric with respect to 1 ^ 2. It vanishes for integer 
rii < OT n2 < 0, because then it becomes the massless vacuum diagram (Fig.EJ, possibly, 
with some polynomial numerator. 
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k 




Figure 4: One-loop massless propagator diagram 




Figure 5: Massless vacuum diagram 



Using Wick rotation and a parametrization ()2.H|1 . we rewrite the definition ()2.19|) of 

G{ni,n2) as 



G{ni,n2) 



IT 



-d/2 



T{m)Tin2) 



(2.20) 



We want to separate a full square in the exponent; to this end, we shift the integration 
momentum: 



k' = k + 



Oil 



-P 



and obtain 

G{ni,n2) 



-d/2 



r(ni)r(n2) 
1 



1^ exp 


aia2 


ai + a2_ 


1^ exp 


a\a2 





n\ — \„ni — \ 



a^ a 



f2.211 



+ a2)-'^/^a1'-^a^^-^dai da2 . 



r(ni)r(n2) 

Now we make the usual substitution ai = rjx, a.2 = ?7(1 — x), and obtain 

1 ' 



G{ni,n2) 



r(ni)r(n2) Jo 
r{-d/2 + ni + n2) 







X 



(2.22) 



r(ni)r(n2) JO 
This integral is the Euler B function, and we arrive at the final result 

r{-d/2 + ni+ n2)T{d/2 - ni)r{d/2 - n2) 



G{ni,n2) 



T{ni)T{n2)T{d - rii - ^2) 



(2.23) 
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For all integer ni^2, these integrals are proportional to 

with coefficients which are rational functions of d. 

The denominator in ()2.19|) behaves as _^ Therefore, the integral 

diverges if c? > 2{ni + 722). At c? — 4 this means ni + n2 < 2. This ultraviolet divergence 
shows itself as a 1/e pole of the ffist F function in the numerator of p.23|l for 711=^2 = 1 
(this F function depends on rii + n2, i.e., on the behaviour of the integrand at — 00). 
The integral ()2.19j) can also have infrared divergences. Its denominator behaves as (A;^)"^ 
at A; — > 0, and the integral diverges in this region if ci < 2?t,2. At d ^ 4 this means n2 > 2. 
This infrared divergence shows itself as a l/e pole of the third F function in the numerator 
of ()2.23|) for n2 > 2 (this F function depends on n2, i.e., on the behaviour of the integrand 
at A; — > 0). Similarly, the infrared divergence at k + p —>■ appears, at ci ^ 4, as a pole of 
the second F function, if ni > 2. 




Figure 6: Complex plane 
Let's consider the integral (Fig. 0} 

as a function of the complex variable p"^. It has a cut along the positive half-axis p^ > 0, 
starting at the branching point = at the threshold of the real pair production (Fig. IH)). 
Analytically continuing it from a negative p"^ = —s (where the function is regular) along 
the upper contour in Fig. IHl p^ = — sexp(— ia) with a varying from to vr, we obtain 

J(s + iO) = Gis-V"". (2.26) 

Similarly, continuing along the lower contour, we get I{s — iO) having the opposite sign of 
the imaginary part, and the discontinuity across the cut is 

J(s + iO) - I{s - iO) = Gxs-m sin(7re) . (2.27) 

At e — > 0, G\ (j2.24j) has a 1/e pole, and the discontinuity is finite: 

/(s + zO) - /(s - iO) = 27rz . (2.28) 
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3 QED at one loop 

3.1 Lagrangian and Feynman rules 

First we shall discuss quantum electrodynamics with massless electron. This theory has 
many similarities with QCD. Its Lagrangian is 

L = ^lJo^lpi^o-\Fo^..F^\ (3.1) 

where ijjQ is the electron field, 

Df,i)Q = {d^ - icqAq^) ipQ (3.2) 
is its covariant derivative, Aofj_ is the photon field, and 

Fq^u = dfj^Aoi, — diyAo^ (3.3) 

is the field strength tensor. 

As explained in the textbooks, it is not possible to obtain the photon propagator from 
this Lagrangian: the matrix which should be inverted has zero determinant. This is due 
to the gauge invariance; in order to have a photon propagator, we have to fix the gauge in 
some way. The most popular way (called the covariant gauge) is to add the term 

= -7^ (d.A^of (3.4) 

to the Lagrangian (j3.ip (oq is the gauge-fixing parameter). This additional term does not 
change physics: physical results are gauge- invariant, and do not depend on the choice of 
the gauge (in particular, on oq). 

= iSo{p) 



P 

p ^ 



Figure 7: QED Feynman rules 
Now we can derive the Feynman rules (Fig. Ej). The free electron propagator is 

^o(p) = i = 4' (3-5) 
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and the free photon propagator is 



DIM = 



(^ \Vi,Vv 
QiJiv - (1 - oo) — 17- 



(3.6) 



The Lagrangian ()3.H) . ()3.4j) contains the bare fields ipQ, Aq and the bare parameters Cq, 
Qq. The renormahzed quantities are related to them by 

ipo = ^y^V" , = Z^l^A , ao = ZAa , cq = ^'^^e (3.7) 

(we shall see why the same renormalization constant Za describes renormalization of both 
the photon field A and the gauge- fixing parameter a in Sect. 13. 2j) . We shall mostly use the 
minimal renormalization scheme (MS). In this scheme, renormalization constants have the 
structure 

Z,(a.) = l + £lii+fff + fHl)fii)V... (3.8) 



£ 47r \ e J \AttJ 

They contain neither finite parts (e^) nor positive powers of e, just terms necessary 
to remove divergences in renormahzed (physical) results. Sometimes, other renormal- 
ization schemes are useful, see Sect. El 

In order to write renormalization constants ()3.8|) . we must define the renormahzed cou- 
pling a in such a way that it is exactly dimensionless at any d. The action is dimensionless, 
because it appears in the exponent in the Feynman path integral. The action is an integral 
of L over d-dimensional space-time; therefore, the dimensionality of the Lagrangian L is 
[L] = d (in mass units). From the kinetic terms in the Lagrangian p.ip we obtain the 
dimensionalities of the fields: [Aq\ = 1 — e, [^/'o] = 3/2 — e. Both terms in the covariant 
derivative ()3.2|1 must have the same dimensionality 1, therefore [eo] = s. In order to define 
a dimensionless coupling a, we have to introduce a parameter with the dimensionality 
of mass (called the renormalization scale): 

<M^,-^^-., (3.9) 

where 7 is the Euler constant. In practice, this equation is more often used in the opposite 
direction: 



(^ = /i^^^^.(«(/^))e-. (3.10) 

We first calculate some physical quantity in terms of the bare charge eo, and then re-express 
it via the renormahzed q;(/x). 

3.2 Photon propagator 

The photon propagator has the structure (Fig. |H)) 

-iD,,{p) = - iDlM + (-z)DjJp)ffl"^(p)(-z)DO,(p) 

+ {-^)DlM^^''^P){-^)Dl,{p)^n'\p){-^)Dl{p) + ■■■ 
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where the photon self-energy ill^^{p) (denoted by a shaded blob in Fig. ^ is the sum of 
all one-particle-irreducible diagrams (diagrams which cannot be cut into two disconnected 
pieces by cutting a single photon line), not including the external photon propagators. 



-)- 



+ 



Figure 8: Photon propagator 

The series fl3.11|) can also be rewritten as an equation 

D,Ap) = DIM + Dl^{p)WP{p)DpM ■ 
In order to solve this equation, let's introduce, for any tensor of the form 



the inverse tensor 

satisfying 
Then 



4-1 - 4-1 



p2 



9 ' 

-i Pf^Pu 
I p2 



D;l{p) = iDXlip)-^M. 
As we shall see in Sect. 13. 3| 

n^.(pK = o, n^,(pK = o. 

Therefore, the photon self-energy has the form 

^f^Ap) = {p'^9t^u - PMPi^)n(p^) , 
and the full photon propagator is 



D^uip) 



p2(i - n(p2)) 



PjiPv 
p2 



Pt^Pu 



(3.12) 



(3.13) 
(3.14) 

(3.15) 
(3.16) 



Its longitudinal part gets no corrections, to all orders of perturbation theory. The full bare 
propagator is related to the renormalized one by 



Therefore, 



DUp;f^)=Dlip';fi) 



p2 



, / ^P^^Py 



(3.17) 
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The minimal ()3.8p renormalization constant ^yi(a) is constructed to make 
finite at e — 0. But the longitudinal part of ()3.17|) containing 

must be finite too. This explains why Za appears twice in ()3.7|) . in renormalization of Aq 
and of qq. 

3.3 Ward identity 

It is easy to check by a direct calculation that (Fig. ^ 



iSoip') ^eo^ iSo{p) = cq [iSo{p') - iSo{p)] 



(3.18) 



Here q = p' — p; substituting ^ and using So{p) = 1/|^, we immediately obtain ()3.18p . 
We shall use graphical notation (Fig. a photon line with a black triangle at the end 
means an external gluon leg (no propagator!) contracted in its polarization index with the 
incoming photon momentum q. A dot near an electron line means that its momentum is 
shifted by q, as compared to the diagram without the longitudinal photon insertion. For 
an infinitesimal q, we obtain from ()3.18p a useful identity 



dSojp) 



-Soiph^,So{p) . 



(3.19) 



p p 



eo 



p 



1r 



Figure 9: Ward identity 
Let's calculate Il^i,{p)p'^ ()3.14|) at one loop: 



eo 



(3.20) 



Two integrals here differ by a shift of the integration momentum, and hence are equal. 

Now we'll do the same at two loops. All two-loop diagrams for Il^^{p)p'-' can be obtained 
from one two-loop diagram with a single photon leg, and using the Ward identity (Fig. E)) 
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we obtain 




(3.21) 



All diagrams cancel pairwisc, except the first one and the last one; these two diagrams 
differ only by a shift in an integration momentum. 

It is clear that this proof works at any order of perturbation theory. 
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3.4 Photon self-energy 

Now we shall explicitly calculate photon self-energy at one loop (Fig. ITO|) . The fermion 
loop gives the factor —1, and 

i{P^9^.u - P^Vumv") = - j Tr le^^^i^^^ie^iu^ • (3.22) 

To simplify finding the scalar function n(p^), we contract in /i and v. In dimensional 

space-time 

5^^ = d, (3.23) 



and we obtain 



^P) ^d-l){-p^)J{2'nY[-{k+pY]{-k^)- ^ 




Figure 10: One- loop photon self-energy 

Now we make a short digression and discuss 7 matrices in d dimensions. Their defining 
property is 

+ 7 V = '^9^" ■ (3.25) 

Therefore, 

1^.1^ = d. (3.26) 
How to find 7^(^7^? We anticommute 7^^ to the left: 

Tm^t'^ = 7m(-7V + 2a^) = -{d - 2)^ . (3.27) 

Similarly, 

^^(/t,l^^' = 7^^(-7''^ + 26^) = {d- 2)(f4 + 21/xIl = Aa ■ b + {d ~ 4:)<lil/> , (3.28) 

and 

^^^^^7/^ = ^^^^[-Y^ + 2c^) = -4a ■ 6 - (rf - 4)^^^ + 2^^^ = -2^^^ - [d - . (3.29) 
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The usual convention is^ 

Trl = 4, (3.30) 

then the usual formulas for traces of 2, 4,. . . 7 matrices hold. 
Now we can find the trace in the numerator of ()3.24|) : 



2, d-2 iel f d'^k 4(fc + p) -k 



^ d-l-p^ J {2TTY[-{k+p)^]{-k^)' ^^'^^^ 

We can set — = 1; it is easy to restore the power of — by dimensionality at the end of 
calculation. The denominators ()2.19|) are 

D^ = -{k + p)\ D2 = -k\ 

All scalar products in the numerator can be expressed via the denominators: 

p' = -l, k^ = -D2, p-k = ^il + D2-D^) . (3.32) 



^rf-2,2 f d'^k -2^2 + 1 + ^2-^1 



We obtain 

n(p^) = 2j-fei J (f^ - — - . (3.33) 

Terms with Di or D2 in the numerator can be omitted, because integrals with a single 
massless denominator (Fig. ^ vanish. 

Restoring the power of —p"^, we arrive at the final result 

n(/)--^2^G.. ,3.34, 

or, recalling (j2.24|) . 

n(/) = ^oi-p2~\ /^-^ . (3.35) 

' (47r)'^/2 (^d-l){d-3){d-4y ^ ' 

3.5 Photon field renormalization 

The transverse part of the photon propagator (j3.16|) is, with the one-loop accuracy. 

Re-expressing it via the renormalized a{fi) (j3.1Uj) . we obtain 

p'D^ip') = 1 + 4^e-^^e^^(7i 4 " ^ 



T 1 ~P 

L = log 



^It would be more natural to use 2'*/^, because 7 matrices can be defined in any even-dimensional space, 
and they are 2''/^ x 2**/^ matrices; moreover, with this definition Fierz identities can be generalized to d 
dimensions. 
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n=2 



We want to expand this at £ ^ 0. Using 

r(l + e) = exp 
where 7 is the Euler constant, and 

00 ^ 

k=l 

is the Riemann ( function, 

C2 = ^, Cs- 1.202, C4 = ^ 
D 90 



(3.38) 



(3.39) 



(3.40) 



we see that 



e-^'g, = 1 + 0{e' 



This is exactly the reason of including exp (—75) into the definition ()3.9|) . 
We obtain 

p D±{p ) = !--■ 



3 47re 



(3.41) 



This should be equal to Z^(a(/x))p^Z}^(p^; /i), where ^^(0) is a minimal ()3.8j) renormal- 
ization constant, and D^(p^; /i) is finite at e ^ 0. It is easy to see that 



ZA(a) 



4 a 



3 Aire 



and 



4a(/i) 



3 Att 



3 



(3.42) 



(3.43) 



The bare propagator D^ij)^) = Z^(a(yu))D^(p^; /i) does not depend on /i. Differentiat- 
ing it in log fi we obtain the renormalization group (RG) equation 



dDl{p^;i2) 



+ ^A{a{fi))Dl{p';fi)=0 



(9 log yU 

where the anomalous dimension is defined by 

d\ogZA{a{^)) 



(3.44) 



7a (a (/i)) 



For any minimal renormalization constant 



dlogfi 



(3.45) 



a 



Zi{a) = 1 + 2:1— + 
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using ()3.9j) 

d\oga{fi) 



-2e + -- 



dlogfi 

we find the corresponding anomalous dimension 

d\og Zi{a{fi)) a{fi) 

7i a /i = = 70— — + ■ • • 

a log /i 4:71 

to be 

7i(a) = -22;i— + ■■■ 
47r 

In other words, the renormalization constant with the one-loop accuracy is 

ZAa) = 1 h ■ ■ • 

The anomalous dimension of the photon field is, from ()3.42|) . 

^Aia) = -- + ■■■ (3.46) 

The dependence of the renormalized propagator (j3.43|) on L can be completely deter- 
mined from the RG equation. We can rewrite it as 

Solving this equation with the initial condition 

pD^{p;f, =-p) = l--—- 

at L = 0, we can reconstruct ()3.43|1 . 

We can also write the RG equation for a(/i). The bare gauge- fixing parameter oq = 
Zyi(a(/i))a(/i) does not depend on fi. Differentiating it in log/x, we obtain 

-^^ + 7A(«(/i)M/i)=0. (3.48) 
a log n 

3.6 Electron propagator 

The electron propagator has the structure (Fig. ITT|l 

iS{p) = z5o(p) +^5o(p)(-OS(p)25o(p) +25o(p)(-z)S(p)25o(p)(-z)S(p)^5o(p) + ■ ■ ■ (3.49) 

where the electron self-energy —iTi{p) is the sum of all one-particle-irreducible diagrams 
(diagrams which cannot be cut into two disconnected pieces by cutting a single electron 
line), not including the external electron propagators. 
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Figure 11: Electron propagator 



The series (I3.49|) can also be rewritten as an equation 

S{p) = So{p) + So{p)J:{p)S{p). (3.50) 

Its solution is 

Sip) = y (3.51) 

Electron self-energy S(p) depends on a single vector p, and can have two 7-matrix struc- 
tures: 1 and 1^. When electron is massless, any diagram for S contains an odd number of 
7 matrices, and the structure 1 cannot appear: 

S(p) =^Sv.(p'). (3.52) 

This is due to helicity conservation. In massless QED, the electrons with helicity A = =f|, 

/ , 

Wl,r = — — V , 

cannot transform into each other. Operators with an odd number of 7 matrices, like ()3.52j) . 
conserve helicity, and those with an even number of 7 matrices flip helicity. Therefore, the 
massless electron propagator has the form 

Sip) = 1 ^ / 2\ 2 ■ (3.53) 

1 - i:vip^)p 




Figure 12: One-loop electron self-energy 
Let's calculate electron self-energy at one loop (Fig. IT^: 
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where we introduced the notation 

^ = l-ao. (3.55) 
To find the scalar function Sy(p^), we take ^ Tr|^ of both sides: 



iel f d'^k N 



(3.56) 



Using the "muhiphcation table" ()3.32|) we obtain 

N = -^ Tr ^7m(I^ + + ^^\ T^^m^ + M 

= -{d-2){p^ +p-k) + ^[ep-k + 2{p-kY -p^P] (3.57) 

where terms with Di or D2 in the numerator were omitted. Restoring the power of — p^, 
we have 

^vip') = - '°(4~)Z'l - 2 - 0^(1, 1) + ^G{1, 2)] . (3.58) 
It is easy to derive a nice property 

G{ni,n2 + 1) {d — 2ni — 2n2){d — ni — ^2 — 1 



G{ni,n2) n2{d-2n2-2 
from the definition (|2.23|) of (^(n 1,722). In particular, 

G(l,2) 



(3.59) 



G(l,l) 

and we arrive at 



(rf-3), (3.60) 



Note that the one-loop electron self-energy vanishes in the Landau gauge Oq = 0. Recall- 
ing ()2.24|) . we can rewrite the result as 

V ( 2^ eli-p')-' d-2 

The electron propagator (j3.53|) . with one- loop accuracy, expressed via the renormalized 
quantities a(/i) (|3.10|) and a(/i) (j3.7p (we may take a(/i) = Oq, because it only appears in 
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the a correction) is 



4vr ^"^(d-3)(c?-4) ^3g3^ 

47re 



It should be equal to Z^{a{fi),a{fi));^Sr{p] fi), where Z^{a) is a minimal ()3.8p renormal- 
ization constant, and Sr{p', /i) is finite at £ ^ 0. It is easy to see that 

Zj,(a, a) = 1 — a (3.64) 

and 

^Sr{p;fi) = l + a{fi){L-l)^. (3.65) 

The bare propagator S{p) = Z^{a{^))Sr{p; jj) does not depend on /i. Differentiating it in 
log n we obtain the RG equation 

^^r(p, /i) ^ ^ a{fi))Sr{p; /i) = , (3.66) 

a log /i 

where the anomalous dimension 

/ / x / NX c;iogZ^(a(/x),a(/i)) 

7^a/i,a/i = — (3.67) 

a log /i 

is ^ 

7^(a, a) = 2a-; h ■ • ■ (3.68) 

47r 

3.7 Vertex and charge renormalization 

Let the sum of all one-particle-irreducible vertex diagrams, i.e. diagrams which cannot be 
cut into disconnected pieces by cutting a single electron or photon line, not including the 
external propagators, be the vertex ieor^{p,p') (Fig. [TSj) . It can be written as 

T''{p,p')=Y + A'{p,p'), (3.69) 

where K^{p,p') starts from one loop. 

When expressed via renormalized quantities, the vertex should be equal to 

r'^ = ZrF^f , (3.70) 

where Z^ is a minimal ()3.8p renormalization constant, and the renormalized vertex F^^ is 
finite at £ ^ 0. 

In order to obtain a physical scattering amplitude (S'-matrix element), one should 
calculate the corresponding vertex (one-particle-irreducible, external propagators are not 
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Figure 13: Electron-photon vertex 
included) and multiply it by the spin wave functions of the external particles and by the 

1/2 

field renormalization constants Z/ for each external particle i. We can understand this 
rule in the following way. In fact, there are no external legs, only propagators. Suppose we 
are studying photon scattering in the laboratory. This photon has been emitted somewhere. 
Even if it was emitted in a far star (Fig. IT^ . there is a photon propagator from the far 
star to the laboratory. The (bare) propagator contains the factor Zj^. We split it into 
Z^"^ ■ ^y^, and put one factor into the emission process in the far star, and the other 

1 /2 

factor Zy( into the scattering process in the laboratory. 




Figure 14: Scattering of a photon emitted in a far star 

The physical matrix element CqTZ^Z^'^ of photon emission (or absorption) by an elec- 
tron must be finite at e ^ 0. Strictly speaking, it is not an S*- matrix element, because at 
least one particle must be off-shell, but this does not matter. This matrix element can be 
rewritten as eTrZa^"^ Z-pZ^Z^'^ . The renormalized charge e and the renormalized vertex 
are finite. Therefore, the minimal ()3.8|) renormalization constant Z]!"^ Z^Z^Z^j^ must be 
finite, too. According to the definition, the only minimal renormalization constant which 
is finite at £ — >■ is 1. Therefore, 

Z^ = {ZpZ^)-''Z-^\ (3.71) 

In order to obtain the charge renormalization constant Z^, one has to find the vertex 
renormalization constant Zr and the electron- and photon-field renormalization constants 
Z^ and Za- In fact, the situation in QED is simpler, because Z-^Z^ = 1, due to the Ward 
identity (Sect. EIS))- 

Starting from each diagram for —iH{p), we can construct a set of diagrams for ieoA'^(p, p'), 
by attaching the external photon line to each electron propagator in turn. Let's calculate 
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the contribution of this set to ieoA^{p,p')q^ using the Ward identity of Fig. El As an exam- 
ple, we consider all vertex diagrams generated from a certain two-loop electron self-energy 
diagram: 



^^^A^^ ^^'Wv^'v^/^^ 



eo 



(3.72) 



eo 



Of course, such cancellations happen for vertex diagrams generated from any self-energy 
diagram (and attaching the external photon line to an electron loop gives 0). Therefore, 
we arrive at the Ward-Takahashi identity 

A'^{p,p')q^ = E{p)-E{p'), (3.73) 

which can also be rewritten as 

T'^{p,p')q, = S-\p')-S~\p). (3.74) 
For g — > 0, we obtain the Ward identity 

A'^(P,P) = -^ or T^ip,p) = ^^. (3.75) 

Rewriting ()3.74|1 via the renormalized quantities, 

ZrT^q, = Z-' [s;\p') - S;\p)] , 

we see that Z^Zr must be finite. But the only minimal renormalization constant finite at 
e ^ is 1: 

Z^Zr = 1 . (3.76) 
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Therefore, charge renormahzation in QED is determined by the photon field renormahza- 
tion: 

Zo. = Z^' . (3.77) 

We know Zr at one loop from the Ward identity ()3.7(i|l and Z^ ()3.(i4j) . Nevertheless, 
let's also find it by a direct calculation. This will be useful, because we'll have to do 
several similar calculations in QCD. We are only interested in the ultraviolet divergence 
of the diagram in Fig. ^1 This divergence is logarithmic. We may nullify all external 
momenta, because terms which depend on these momenta are convergent: 

Of course, we should introduce some infrared regularization, otherwise this diagram van- 
ishes. We have 



Averaging over k directions: 
we obtain (4-dimensional 7-matrix algebra may be used) 



A" = --^ao7"/ (0,(Z^- (3.80) 




Figure 15: One- loop QED vertex 

Now let's find the ultraviolet divergence (1/^) of this integral, introducing a sharp 
infrared cutoff into the Euclidean integral (we may use ()2.12|1 here) 



d'^k 



{2TTY{-k^y 



uv 



87r2 



k 



dk 



-2e 



(47r)% (47r)2£ 



(3.81) 



Any infrared regularization can be used; instead of a cut-off, we could insert a non-zero 
mass, for example: 



d'^k 



(27r)^(-P)^ 



uv 



d'^k 



im 



-2e 



{2ttY (m2 - P)2 (An)' 



-Tie] 



(47r)2e 



(3.82) 
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(see dUD, (EID). 

Therefore, the Xje part of A" is 

A" = (3.83) 

(it does not depend on momenta). This calculation can be viewed as the simplest case of 
the infrared rearrangement [2j. We obtain the renormalization constant 

Zr = 1 + a-^ , (3.84) 
47re 

in agreement with ()3.76p . ()3.64|) . 

1/2 jT^^rm 

The bare charge Cq = Za {a{^))e{fi) does not depend on /i. Differentiating ()3.10p in 
dlogfi, we obtain the RG equation for a{fi): 

(iloga(u) ^/ / XX /X 

= -2e - 2/3 a /i , 3.85 

where /3 function is defined by 

Of ( 1 ^ log («,(/!)) 

/3 a, = — — . 3.86 

/ a log jjL 

With one-loop accuracy, differentiating 

Za{a) = 1 + zi- \ 

47r£ 

we can retain only the leading term —2e in ()3.85|) : 

P{a)= (3o— + ■■■ = -z^— + ■■■ 
A^[ An 



In other words, at one loop 



a 

Z„(a) = l-/3o— + ••• (3.87) 



We have obtained (see (jTZ?)), (HO^ 

ZT^ = 1 + 

3 47re 



Z^ = Zf = l + \^ + --- (3.88) 



Therefore, the QED (3 function is 

/3(a)=/3o£ + --- = -^£ + --- (3.89) 

When we consider the renormalized a;(/i) after taking the physical limit e — > 0, we may 
omit -2e in dUSl): 

^^^ = -2/J(a.W). (3.90) 
a log jj, 
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In QED Po = —4/3, and the /3-function ()3.89p is negative, at least at sufficiently small 
a, where perturbation theory is vahd. Therefore, the running grows with fx. This 
corresponds to the physical picture of charge screening: at larger distances (smaller /i) the 
charge becomes smaller. 

The RG equation ()3.9()|1 in the one-loop approximation. 



can be rewritten as 



and its solution is 



dlogfj, 4:71 \ An 

d An 



dlogjjL a{jj) 



2(3o , (3.92) 



An An ^ , a' 



Finally, 



a(/i') a(/i) /i 

"^^'^ = ^) — 7 ■ ^^-^^^ 

l + 2/3o^log^ 
An jji 



3.8 Electron mass 



Until now we treated electrons as massless. This is a good approximation if characteristic 
energies are large (in the Z region at LEP, say). Now let's recall that they have mass. The 
QED Lagrangian ()3.1|) now is 

L = V-o {lip - mo) ^0 + ■ ■ • (3.94) 
and the free electron propagator (|3.5|) becomes 

So{p) = -1 = — 2 • 3.95 

p — rriQ — mo 

In addition to ()3.7|) . we have also mass renormalization 

mo = Zm{a{ii))m{ii) . (3.96) 

The electron self-energy has two 7-matrix structures, because helicity is no longer con- 
served: 

S(p)=|^Sv(p')+moS5(/), (3.97) 

and the electron propagator is 

^^^^ ^ i> ~ mo - i^llvip') - mo-Lsip') ^ 1-Sy(p2) l + T.s{p^) ~ ' ^^'^^^ 

1 - Sy(r) 
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It should be equal to Z^Sr{p;fi)- The renormalization constants are determined by the 
conditions 

1 + S 

(1 — Ey)Z^ = finite , :^Zm = finite . (3.99) 

1 — Ly 

The first of them is the same as in the massless case; the second one, defining Z^, can be 
rewritten as 

(1 + J:s)ZfZm = finite . (3.100) 
The one-loop electron self-energy (Fig. IT^ is 

In order to single out S5, we should retain the mo term in the numerator of the electron 
propagator (it fiips helicity). At large p"^, 



= -iel{d-Cl 



{2'kY {k + pfk^ V k^ 
d'^k 1 



{2^Y{k + pYk'^ ' 

Retaining only the ultraviolet divergence ()3.8H) and re-expressing via renormalized quan- 
tities (this is trivial at this order), we obtain 

Ss = (3 + a(^))^. (3.103) 
From ()3.100p we see that a terms cancel: 

(1 + T.s)Z^Zm = (1 + (3 + a)£^) (1 - a^) = 1 , 
and Zm is gauge-independent: 

Z^ = l-3-^ + --- (3.104) 

A-KE 



The bare mass mg = Zm{oi{jji))m{n) does not depend on /i: 

dm{n) 



d\og /X 

where the mass anomalous dimension is 



+ 7™(a(/i))m(/i) = 0, (3.105) 



From (j3.104j) we obtain 



rflogZ„(a(/i)) 

7m(a(/^)) = -T. • (3.106) 

d log jjL 



a a 

lm[a) = Imo-r = ^ir~ ~^ (3.107) 

47r 47r 
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In order to solve the RG equation ()3.105p for let's write it down together with the 

RG equation ()3.85p for a{fi): 

d^oga dlogm 

— = -2/3(a) , — = --fm{a) , 

a log /i a log n 

and divide the second equation by the first one: 

dlogm _ 7n^(a) 
d\oga 2/5(a) 

The solution is 

/■"('"'^ ^ (a) da 

m(/i') = exp / (3.109) 



(3.108) 



At one loop, we obtain from ()3.89|) and ()3.107|) 

m{fi) =m{fi) I—- \ =m(/i) — - . (3.110) 

Va(/i) J J 

The running electron mass m(yu) decreases with /i. 

4 QCD at one loop 

4.1 Lagrangian and Feynman rules 

The Lagrangian of QCD with nf massless fiavours is 

L = J2<lo^^Jpqo^-\G^o,.G^, (4.1) 

i 

where goi are the quark fields, 

D^.qo = (^M - ^^?o^om) % , = A^o^t'' (4.2) 

are their covariant derivatives, A^^ is the gluon field, are the generators of the colour 
group, and the field strength tensor is defined by 

[D^, D^]qo = -igoGo^^uqo , Gq^^ = G^^^X ■ (4-3) 

It is given by 

G^o,. = d,Al - d^Al^ + g.r^'Al^Al^ , (4.4) 
where the structure constants j""^" are defined by 

\tX\=ij^^H\ (4.5) 
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Due to the gauge invariance, it is not possible to obtain the gluon propagator from this 
Lagrangian. We should fix the gauge. The most popular way (called the covariant gauge) 
requires to add the gauge-fixing term (with the parameter ao) and the ghost term to the 
Lagrangian ()4.H1 : 

AL = {d.Al^f + {d^cl){D,cl) . (4.6) 

Here Cq is the ghost field — a scalar field obeying Fermi statistics, with the colour index a 
(like the gluon). Its covariant derivative is 

D,cl = {8,6'^' - z^oO 4 , K = "^o.itT , (4.7) 

where 

are the generators of the colour group in the adjoint representation. 

= iiS'o(p) 



V 

fx p V fiiyyi/ 



P 



Figure 16: Propagators in QCD 

The quadratic part of the QCD Lagrangian gives the propagators shown in Fig.^J The 
quark propagator is So{p) ()3.5p : the unit colour matrix is assumed. The gluon propagator 
is D°^(p) (j3.6p : the unit matrix in the adjoint representation 6'^'' is written down explicitly. 
The ghost propagator is 

Goip) = \ (4.9) 



times 5"^. 

/i a 



t°- X igo'j^ 



Figure 17: Quark-gluon vertex 

The quark-gluon vertex is shown in Fig. El The problem of calculation of a QCD 
Feynman diagram reduces to two separate subproblems: 

• calculation of the colour factor; 
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• calculation of the "colourless" diagram. 

The first subproblem can be formulated as calculation of the colour diagram. It looks like 
the original QCD diagram, where quark lines and gluon lines mean the unit matrices in the 
fundamental representation and in the adjoint one; ghost lines can be replaced by gluon 
lines in colour diagrams. The quark-gluon vertex in a colour diagram gives t*^, the first 
factor in Fig. El When calculating the "colourless" diagram, all factors except the colour 
structure are included; the quark-gluon vertex gives igo'j^, the second factor in Fig. IT71 




Figure 18: 3-gluon vertex 

The 3-gluon vertex is shown in Fig. ^1 When separating it into the colour structure 
(which appears in the colour diagram) and the Lorentz structure (which appears in the 
"colourless" diagram), we have to choose some rotation direction. These rotation directions 
must be the same in the colour diagram and in the "colourless" one. The colour structure 
is if^-i^-^a-i^ where the colour indices are written in the chosen order (in Fig. clockwise). 
The Lorentz structure is 

^mM2M3(p^^p2^p3) = (p3 - p2)^^^7^^^^ + (pi - P^Tg^'^' + (P2 - PiY'g^'^' , (4.10) 

where the polarization indices and the momenta are numbered in the same way, here 
clockwise. Of course, the full 3-gluon vertex does not know about the rotation direction: 
if we reverse it, both the colour factor and the V tensor change sign, and their product 
does not change. 




Figure 19: 4-gluon vertex 

The 4-gluon vertex contains terms with 3 different colour structures. Therefore, any 
diagram containing at least one 4-gluon vertex cannot be factorized into the colour diagram 
and the "colourless" one. This is very inconvenient for programs automating diagram 
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a 



II, rv 




Figure 20: Propagator of the auxiliary field 



a 



b 




Figure 21: Vertex of the auxiliary field interaction with gluons 



calculations. The authors of (at least) two such programs (CompHEP and GEFICOM) 
invented the same trick to ensure the possibility of calculating the colour factor as a 
separate subproblem. Let's say that there is no 4-gluon vertex in QCD, but there is a 
new field interacting with gluons (this field is shown as a double zigzag line in Fig. IT^ . 
This is an antisymmetric tensor field; its propagator is shown in Fig. 1201 and the vertex of 
its interaction with gluons — in Fig. |^ Its propagator in the momentum space does not 
depend on p. Therefore, in the coordinate space, it contains S{x) — this particle does not 
propagate, and two vertices of its interaction with gluons (Fig. IT!I|) are at the same point. 
In the colour diagram, lines of this particle can be replaces by gluon lines, and vertices — 
by 3-gluon vertices. In "colourless" diagrams, the second factor in Fig. |^ is used for the 
vertices. The rotation directions must agree, as usual. 



Finally, the ghost-gluon vertex is shown in Fig. |221 It contains only the outgoing ghost 
momentum, but not the incoming ghost momentum. This is because its Lagrangian ()4.fi|l 
contains the covariant derivative of c but the ordinary derivative of c. In colour diagrams, 
we replace ghost lines by gluon lines, and vertices — by 3-gluon vertices; their rotation 
direction is fixed: the incoming ghost the outgoing ghost the gluon (Fig. |221)- 

Calculation of colour diagrams is discussed in Appendix 1X1 

The renormalized fields and parameters are related to the bare ones, similarly to 



Qio = Zg/'^qi , Ao = Z^aA , Co = Z\l'^c , ao = ^ao , 5^0 = Z]l'^g ; (4.11) 



Ai b 




V 



Figure 22 



Ghost-gluon vertex 



QED dSID, by 
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the QCD running coupling as{^) is 
4.2 Quark propagator 

The quark propagator in QCD has the same structure ()3.49|) as the electron propagator in 
QED (Sect. 13. 6p . In massless QCD, the quark self-energy has the form 

Eip) = ^Evip') . (4.13) 

The one-loop diagram (Fig. 0^ differs from the QED one ()3.6ip only by the replacement 
Co go and by the colour factor Cp ()A.23|) : 

S^v(P^) = -C.^^|^^a,G, . (4.14) 

Therefore, the quark field renormahzation constant with one-loop accuracy is 



Zg = l-CFa-^ + --- (4.15) 
Aire 



and the quark-field anomalous dimension is 



^ 2CFa-^ + --- (4.16) 
47r 




Figure 23: One-loop quark self-energy 

If the quark flavour gj has a non-zero mass, its self-energy has also the helicity-fiipping 
structure ^^(p^). The calculation of its one-loop ultraviolet-divergent part (Sect. 13. 8|) 
applies in QCD practically unchanged, with the substitution a — > and the extra colour 
factor Cp'- 

S5 = C^(3 + a)-^ + --- (4.17) 

Ane 

Therefore, the quark-mass renormahzation constant is 

Z^ = 1-?,Cf-^ + --- (4.18) 
47re 



and the mass anomalous dimension is 



ry 

lm = QCF-^ + --- (4.19) 

47r 
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4.3 Gluon propagator 

Ward identities in QCD are more complicated than in QED. Nevertheless, the gluon self- 
energy iS°'''Ilij,i^{p) is transverse: 

n^.(p)p^ = , n^,(p) = {p^g^, - p^p.Mp^) , (4.20) 

as in QED. Explanation of this is beyond the level of these lectures. 




Figure 24: One-loop gluon self-energy 

At one loop, the gluon self-energy is given by 3 diagrams (Fig. The quark-loop 
contribution differs from the QED result ()3.34|1 only by the substitution eo ^ go, and the 
colour factor TpUj: 




(4.21) 



Figure 25: Gluon-loop contribution 

The gluon-loop contribution (Fig. has the symmetry factor ^ and the colour factor 
Ca ()A.29|) . In the Feynman gauge ao = 1, 



2 {2tiY k^{k + pf ' (4.22) 

= V^^p{p, -k - p, A;)V^^^"(-p, -k, k+p). 

The Lorentz part of the 3-gluon vertex (j4.1Up here is 

Vf^apip, -k -p,k) = {2k + p)^gap - {k - p)agpfi -{k + 2p)pg^a ■ 
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V^^^°'(— p, —k, k + p) coincides with it. Therefore, the numerator in ()4.22|) is 



N = d [{2k + pf + {k~ pf + {k + 2pf] 

- 2{2k + p)-{k-p)- 2{2k + p) ■ {k + 2p) + 2{k + 2p) ■ {k - p) . 

Using the "multiphcation table" ()H.H2j) and omitting terms with Di^2 in the numerator 
(which give vanishing integrals), we have 



p^ = -l^ e = -D2^Q, p-k = -{l~Di + D2 



Finally, we arrive at 



1 

2 ' 



(4.23) 



k 



"k^p 



Figure 26: Ghost-loop contribution 

The ghost-loop contribution (Fig. has the fermion-loop factor —1 and the colour 
factor Ca <\EM- 



iCAgo 



d'^k k-{k + p) 1 Qoi-p^) 



2\l-e 



{27TYk^{k+py 



(4vr)rf/2 



-Gi. 



(4.24) 



These two contributions, taken separately, are not transverse. Their sum is transverse, 
and hence has the structure ()4.2()j) : 



(rf-l)(-p2) 

In an arbitrary covariant gauge. 



(477)^/2 '2{d-\)' 



(4.25) 



n,(p^) = Ca^ 
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{AttY/^ 2{d-l) 



3d-2 + {d- l){2d - 7)e - -{d -l){d- 4)e^ 



(4.26) 
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(here ^ = 1 — ao fj3.55|) : we leave the derivation as an exercise for the reader). 

The transverse part of the gluon propagator with one-loop accuracy, expressed via 
renormalized quantities and expanded in e, is 



1 + ^Me-L- 



1/ 13\ 4 

a ^ I Ca -LpUf 



+ 1 CA--T,n,]e 



where L = log(— Therefore, the gluon field renormalization constant is 

'1 



(4.27) 



1 - 



as 
Atts 



+ ■ ■ 



(4.28) 



We don't write down the expression for the renormalized gluon propagator; it is easy to 
do this. The gluon field anomalous dimension is 



7a = 

Making the substitutions as 

suits dm, (Hini). 



13 



3 y '^A + -TpUf 



4:71 



+ 



(4.29) 



a, Ca — > 0, TpUf — 1, we reproduce the QED re- 



4.4 Ghost propagator 

The ghost propagator has the structure 

iG{p) = iGoip) + tGo{p){-t)i:iGo{p) + iGo{p){-i)^tGo{p){-i)^tGo{p) + ■■■ (4.30) 

where the ghost self-energy — is a scalar function of p^, and the free ghost propagator 
Go{p) is given by ()4.9|1 . Therefore, 



G{p) 



1 



p2 — S(p2 



At one loop (Fig.|2Z|), 

S(P^) 



-iCAgo 



, I d^k p'^jk + pY / k^k, 
{2'kY k\k + pf V^'^" ^ F 



2\l-£ 



Ga 



(47r)"!/2 



-^G(1,1) + ^G(1,2) 



1^ g'oi-p'y-' 

4^ (47r)'^/2 



Gi [d 



{d — 3)ao] . 



(4.31) 



(4.32) 
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Figure 27: One-loop ghost self-energy 



The propagator, expressed via renormalized quantities and expanded in e, is 



G{P) 



^ + ^^1^" 4 



Therefore, the ghost field renormalization constant is 

3 — a a. 



4 Ane 



+ 



and the anomalous dimension is 



7c 



^ 3 — a as 



(4.33) 



(4.34) 



(4.35) 



4.5 Quark— gluon vertex 

Let's find the ultraviolet divergence of the quark-gluon vertex at one loop fFig. OHj). The 
first diagram differs from the QED one only by the colour factor Cp — Ca/'^ (jA.3Hl . and 



A? 



Ca\ CXs 

4:7ie 



(4.36) 




Figure 28: Quark-gluon vertex at one loop 

The second diagram fFig. I29|) has the colour factor C^/2 ()A.33|) . It has a logarithmic 
ultraviolet divergence. Therefore, as in Sect. 13. 7^ in order to find the 1/e term we may 
nullify all the external momenta and provide some infrared cutoff. In the Feynman gauge 
(oo = 1), 



36 



a 




Figure 29: The non-abelian diagram 
where the Lorentz part of the 3-gluon vertex ()4.1U|) is 

Therefore, 



Averaging over k directions: 



tt" ^ —7° 
a 

and using the uhraviolet divergence ()3.81|) . we obtain 



In the arbitrary covariant gauge 

K = \{l + a)CA^r (4.39) 

(derive this result!). 

The 1/e term (the ultraviolet divergence) of the one- loop quark-gluon vertex is, 
from (ICTl) and 

A" = f Cpa + Ca'^\ :^7" • (4.40) 



4 / 47re 

Therefore, the quark-gluon vertex renormalization constant is 

a + 3\ a 



Zr = l+[Cpa + CA^\^^ + --- (4.41) 
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4.6 Coupling constant renormalization 

The coupling renormalization constant is 

Z^ = {ZrZ,)-^Z^\ (4.42) 
similarly to Sect. Kl 71 The product Z^Zq, with one-loop accuracy, is, from ()4.4H] and ()4.15j) . 

ZrZ, = l + c/-±^^ + --- (4.43) 
4 Ane 

In QED, it was equal to 1 ()3.7(j|l . due to the Ward identity (Sect. 13. 7j) . Ward identities 
in QCD are more complicated. Making the replacements as — > a, Cp — 1, Ca 0, we 
reproduce the simple QED result. 

The coupling renormalization constant, with one-loop accuracy, is, from ()4.43|) and ()4.28p . 

It is gauge-invariant. This is a strong check of our calculations^. The /3-function fj3.86|) in 
QCD is, with one-loop accuracy, 

a 114 
/3(a.)=/3o^ + --- (3o = -CA--Tpnf. (4.45) 

If there are not too many quark flavours (namely, for Nc = 3, if nj < 16), then /3o > 0, in 
contrast to the QED case f3o = —4/3. 
The RG equation 

rfloga,(/i) 

— — = -2/? a,^ 4.46 

dlog/i 

shows that as(/u) decreases when /i increases. This behaviour (opposite to screening) is 
called asymptotic freedom. Keeping only the leading (one-loop) term in the /3-function, we 
can rewrite the RG equation in the form ()3.92|) . The running coupling at the renormaliza- 
tion scale /i' is related to that at /i by ()3.93p . but now with positive f3o- The solution of 
the equation ()3.92j) can also be written as 

(^sif^) = (4-47) 

Po log -r^ 

where plays the role of the integration constant. It sets the characteristic energy scale 
of QCD; the coupling becomes small at 3> Aj^. 

Coupling-constant renormalization in the non-abelian gauge theory was first considered 
in [Hj, and the sign of the /3-function corresponding to asymptotic freedom was obtained. 
The magnitude was not quite right (it is right for a spontaneously-broken gauge theory 

•^We have actually presented detailed calculations of Za and Zr only in the Feynman gauge uq — 1; 
this is enough for obtaining 14.44|l . 
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with a higgs). The first correct calculation of the one-loop /3-function in the non-abelian 
gauge theory has been published in [Uj. It was done in the Coulomb gauge, which is 
ghost-free. The contribution of the transverse-gluon loop has the same sign as that of 
the quark loop, i.e., leads to screening. However, there is another contribution, that of 
the loop with an instantaneous Coulomb gluon. It has the opposite sign, and outweights 
the first contribution. This /3-function has also been calculated by 't Hooft but not 
published (mentioned after Symanzik's talk at a meeting in Marseilles in 1972). Later 
it was calculated in the famous papers [TTl [T2]. The authors of these papers applied 
asymptotic freedom to explain the observed behaviour of deep inelastic electron-proton 
scattering. This was the real beginning of QCD as a theory of strong interactions. The 
authors of [IH IT^ received the Nobel prize in 2004. 

4.7 Ghost— gluon vertex 

All QCD vertices contain just one coupling constant. Its renormalization can be found 
from renormalization of any vertex: quark-gluon fSect. 14. 6|) . ghost-gluon, 3-gluon, or 
4-gluon. Here we shall derive it again by calculating the ghost-gluon vertex at one loop. 



The first diagram (Fig. EOl) has the colour factor Ca/2 (IX35|i . The ultraviolet 1/e 
divergence of this diagram is proportional to the bare vertex (i.e., to the outgoing ghost 
momentum p°'); the coefficient diverges logarithmically. We may nullify the external mo- 
menta except p, and use ()3.8H) : 



a 




k 



Figure 30: One-loop ghost-gluon vertex (diagram 1) 



A' 




/ 



{27rY (A;2)3 



d'^k k'^pi'k 



.V 




) 




(4.48) 



Ca 2 1 , 



I 



d'^k 1 1 



{27iY (P)2 8 



Aire 



■p 
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a 




Figure 31: One-loop ghost-gluon vertex (diagram 2) 
The second diagram (Fig. IHTj) has the same colour factor; in the Feynman gauge (oq 



1), 



'^'^ J (2^0^ m ^ 

3 .Ca 2 a f d'^k 1 3^ a 



In the arbitrary covariant gauge, 

= s^^^i;;!^" (4.50) 

(derive this result!). 

The full result for the 1/e term in the ghost-gluon vertex at one loop is 

1 a 

A- = -CAa-^p'' . (4.51) 

2 Ane 

Therefore, the vertex renormalization factor is 

1 a 

Zre = l + -C^a-^ + --- (4.52) 

2 Aire 

The coupling renormalization constant 

= {ZrcZcY^Z^^ (4.53) 

is, from 

ZvcZc = 1 + Ca- 



4 A-KS 

(which follows from (j4.34|) and (j4.52|) ). equal to (j4.44|) . Thus we have re-derived [3q (j4.45j) : 
this derivation is, probably, slightly easier than that from the quark propagator and quark- 
gluon vertex. 
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5 Two-loop corrections in QED and QCD 

5.1 Massless propagator diagram 

Let's consider the two-loop massless propagator diagram (Fig. l32|) . 



-k\-p^Y ^"'G(ni, 712,^3, 714,715) 



(5.1) 



Do 



-kl , = -kl , D, 



The power of —p^ is evident from dimensionality. Our aim is to calculate the dimensionless 
function G{ni,n2, n^, 714, 715). It is symmetric with respect to the interchanges (1 2, 3 <-> 
4) and (1 ^ 3,2 ^ 4). It vanishes when indices of two adjacent lines are non-positive 
integers, because then it contains a no-scale subdiagram. 




ki+p 



k2+p 



Figure 32: Two-loop massless propagator diagram 



n3 



714 




71i 



n2 



Figure 33: Trivial case 715 = 



When one of the indices is zero, the problem becomes trivial. If 715 = 0, it is the product 
of two one- loop diagrams (Fig. ESI): 

G(71i, 7l2, 7I3, 714, 0) = ^(711, 713)^(712, 7I4) . (5.2) 

For integer 7ij, the result is proportional to Gl = G(l, 1, 1, 1, 0), see ()2.24|) . 
If ni = (Fig. I34|l. then the inner loop gives 

G(7l3,7l5) 



.^2-)n3+n.5-rf/2 
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n2 



n4 + ris + ns - d/2 




X 




77,5 



Figure 34: Trivial case rii = 



n2 



and hence 

G(0, n2, n4, ris) = G{n3, n5)G{n2, n^ + us + n^- d/2) . (5.3) 

The cases ^2 = 0, = 0, n4 = are given by the symmetric formulae. For integer rij, the 
result is proportional to 

4^2 



G2 = G(0, 1, 1, 0, 1) = ^(1, 1)G(1, 2 - d/2) 



{d -3)id- 4)(3d - 8){3d - 10) 
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(Fig. ESI). 



T{l + 2e)T\l-e) 
r(l -3£) 



(5.4) 




Figure 35: Basis integral 

But what can we do if all 5 indices are positive? We shall use integration by parts [H] 
— a powerful method based on the simple observation that integrals of full derivatives 
vanish. When applied to the integrand of 15.11 the derivative d/dk2 acts as 

^^^2ik2 + p) + ^2k2 + ^2ik2-h). 

Ok2 D2 -D4 -D5 

Applying {d/dk2) ■ {k2 — ki) to this integrand and using 

{k2 - k,)^ = -D, , 2k2 ■ ik2 - k,) = D,-D^-D,, 
2{k2+p)-ik2-k^)=D,-D2-D,, 



we see that this operation is equivalent to inserting 

n2 



d — n2 — — 2nz 



?T,4 
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under the integral sign in ()5.H) (the term d comes from differentiating k2). The resulting 
integral vanishes. 

On the other hand, we can express it via G with shifted indices. Let's introduce the 
notation 

l^G(ni, n2, ns, ^4, n^) = G{ni ± 1, n2, ris, n^, n^) , (5.5) 

and similarly for 2^, etc. Then the relation we have derived (it is called the triangle 
relation) takes the form 



[cf - ^2 - 714 - 2715 + ^22+(l" - 5") + n44+(3~ - 5")] G = . 



(5.6) 




Figure 36: Integration by parts 

From ()5.6|) we obtain 

^^n.2^(5--l-) + .M^(5--3-)^^ (5.7) 

a — n2 — — 2n^ 

Each application of this relation reduces ni +113 + by 1 fFig. IHUj). Therefore, sooner or 
later one of the indices rii, n^, will vanish, and we'll get the trivial case ()5.2|) . ()5.3|) . or 
symmetric to it. 

Let's summarize our achievements. There is one generic topology of two-loop propa- 
gator diagrams (Fig. IHTj) . This means that all other topologies are reduced cases of this 
one, where some line (or lines) shrink(s) to a point, i.e., the corresponding index (-ices) 
vanish^. All Feynman integrals of this class, with any integer indices ni, ?t,2, 713, 714, 775, can 

''At first sight it seems that the diagram with a self-energy insertion into a propagator (Fig. I38|l is not 
a reduced case of Fig.|^ But this is not so: we can collect two identical denominators together, and then 
it obviously becomes Fig. 1371 with one line shrunk. 
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Figure 37: Two-loop propagator diagram 




Figure 38: Insertion into a propagator 




Figure 39: Basis integrals 



Figure 40: Sunset diagrams 
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be expressed, by a simple algorithm, as linear combinations of two basis integrals (Fig. 
with coefficients being rational functions of d. The basis integrals are G\ and G2-, where 



G 



rd + ne)r"+i('i -£) 



Til - {n + l)e) 



is the n-loop sunset diagram (Fig. HO)). 

Let's consider an example: Fig. EH with all rij = 1, i.e., G(l, 1, 1, 1, 1). Applying ()5.7|) 
to it, we obtain 





(5.9) 



where the symmetry has been taken into account. Here a dot on a line means that this 
denominator is squared. We have 



^(1, 1,1, 1,1) 



d-A 
2 

d-A 
2 

d-A 



[G(l,2,l,l,0)-G'(0,2,l,l,l)] 
[G(l,l)G(2,l)-G(l,l)G(2,3-rf/2)] 



(5.10) 



G(2,l) G(2,-d/2 + 3) ^ 



G(l,l) 



G(l,-d/2 + 2) 



because 



Using (ITHni) . 



Gi = ^(1, 1) , G2 = ^(0, 1, 1, 0, 1) = ^(1, 1)^(1, -d/2 + 2) 



^(2,3-^/2) (3rf-8)(3(i- 10) 



and we obtain 



G(l,l,l,l,l) 



^(1,2-^/2) 



-(^-3)^? + 



d-A 



(3rf-8)(3c/- 10) 



d-A 

It is easy to expand this result in e using 

r(l + 2£)r2(l-2e) 



d-A 



Go 



i2 - gl) 



{d-?,){d-Af 



(5.11) 



92 



g\ r2(i + 5)r(i-5)r(i-3e) 



1 - eCs^' + 



(5.12) 
(5.13) 



We arrive at 

^(1, 1,1,1,1) = 6C3 + --- 

Three-loop massless propagator diagrams also can be calculated using integration by 
parts "B]. 
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5.2 Photon self-energy 



The photon self-energy at two loops is given by 3 diagrams (Fig. l^T|) . It is gauge- invariant, 
because there are no off-shell charged external particles. Therefore, we may use any gauge; 
the calculation in the Feynman gauge ao = 1 is easiest. 




Figure 41: Two-loop photon self-energy 




Figure 42: The first contribution 

The first two diagrams contribute equally; we only need to calculate one (Fig. and 
double the result. This diagram contains one-loop electron self-energy subdiagram —i'^Y^v'- 

Substituting the result 1)3.611) in the Feynman gauge Oq = 1, we obtain 

^ el ^ d-2 f d'^k TT^,ilt + ^)rit 



(47r)'^/2^''^ J {27iY DiDl+' ■ 
Using the "multiplication table" ()3.32|) and omitting Di in the numerator, we have 

ni;: = - ""^fj^^^ d - 2fG{l, 1) [G(l, l+e)- G(l, e)] 

Using the property ()3.59p . we finally arrive at 

e^(_V)^ 2{d-2 ) 
(47r)^ ' rf-4 
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Figure 43: The second contribution 



The second contribution (Fig. is more difficult. It is a truly two-loop diagram 



(Fig. 



(2^(2^ 



Tr ieo'y^.i 



—I 



kf " kl{k,-k,f 



or 



-fin 



cfiki d'^k2 



N 



V i2nri2nrD^D2D,D,D,' 

N = Ti^^ih + i>)iAh + ^h^hYh ■ 

According to (jS.ip . the "multiplication table" of the momenta is 
/ = kl = -Ds, kl = -D,, 

p■k^ = ^il + D■s-D^), p . k2 = ^{1 + D,, - D^) 

ki-k2 = ^{D5-Ds-D,). 



(5.15) 



All products DiDj in the numerator can be omitted (produce vanishing integrals), except 
D1D4 and -D2-D3. Calculating the trace with this simplification, we obtain 



N^2{d- 2) -{d - A){DiDi + D2D3) 

+ 2{D, + D2 + D3 + D4) - 2Dl + {d- 8)^5 - 2 



Therefore, 



(47r)' 



2{d - 2) -2(d - 4)G(0, 1, 1, 0, 1) + 8G{0, 1, 1, 1, 1) 



2G(1, 1, 1, 1, -1) + id- 8)G(1, 1, 1, 1, 0) - 2G(1, 1, 1, 1, 1) 
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Some integrals here are trivial: 

G'(1,1,1,1,0) = G?, ^(0, 1,1, 0,1) 

or simple: 



G(0, 1, 1, 1, 1) = G(l, 1)G(1, 3 - d/2) = ^3jG2 

(here we used fj3.59p ). The only truly two-loop integral, G(l,l,l,l,l), has been calculated 
in Sect. 15.11 and is given by (|5.11|) . But what is (7(1,1,1,1,-1)? This integral would 
factor into two one-loop ones (Fig. IH^ if it had no numerator. In principle, it is pos- 
sible to calculate such integrals using integration-by-parts recurrence relations. But we 
have not discussed the necessary methods. Therefore, we shall calculate this integral in a 
straightforward manner: 



1 r d'^ 

G(U.l.l.-l) = -^/j^ 



[-ih 



]=-- 



ki d'^ki 



k2 d'^k2 



Both vector integrals are directed along p; therefore, we may project them onto p: 

2 f f ki-pd%V 1 



G(l, 1,1, 1,-1; 



'jl-dp2 



D1D3 



(5.16) 



Collecting all this together, we get 



d-A 



(5.17) 



Separate contributions to li^uip) are not transverse, but their sum is, i.e., it has the 
structure ()3.15|1 . The full two-loop photon self-energy 



IS 



et{-p')-'\ 



d 



{AtiY (rf-l)(rf-4) 
The one-loop result (Sect. 13. 4|) is given by ()3.34p . 



2ni;: + Ti2^^ 

(rf-l)(-p2) 



d-A 



(5.1^ 



5.3 Photon field renormalization 

The transverse part of the photon propagator, up to e^, has the form 



p'D^ip' 



1 -i + it^/^(,) + ^o(V)-^^ 



n(p2 



{Aixyi^ 
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(47r) 



f2{e) + 



(5.19) 



From (poll and ^TT^ . using ^TI^ and we have 

l-e 



efiie) = - 



(1 -2£)(3-2e) 



l-e 



^1 



£(3 - 2e) 



6-3e + 4e2-4£3 2 2 - 2e + 
(1 - 25)2(3 - 2e) " ( 



l-3£)(2-3e) 
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We can expand these functions in e: 



ee^'fi{e) = do + cue + cue"^ + 



At one loop, using 



we obtain 



91 = 1- -C2e' + 



4 20 

e + {-C 

3 9 ' 



112 
~27 



At two loops, using ()5.12|) . we see that 1/e terms cancel, and 



2 o-ve ^ / N 16 106 

e f2ie) = \ e 

J2K J g 27 



I6C3 



16 



7 



<2-- e^ + - 



In order to re-express it via a{fi), it is sufficient to use ()3.10|) with the one-loop renor- 
malization constant ()3.87j) . ()3.88|1 : 



(5.20) 



where 



L = I06 



(indicating the argument of a in the term would be beyond our accuracy). Substituting 



the e expansions, we obtain at L = (i.e. /i^ 

p^D^ip^) = 1 + ^^(cio + Cue + cue^ + ■■■) 
Aire 

/ (y \ 2 

+ [-^J [c2o + C2ie + C22£^ H /5o(cio + Cue + 0^2^^ ^ )] + 

This should be equal to Za(«(/x))p^£^^(p^; /i), where 



ZaM = 1 + —zi + 
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(2:20 + 221^) + 



(5.21) 



and 



a 



ra + ■ ■ ■ + 



1 + 

Equating a terms, we obtain 

Zl = Cio , Ti = Cii , rii = Ci2 . 

Equating terms, we obtain 

2:20 = C20 - PoCio , 2:21 = C21 - (cio + Po)cii , r2 = C22 - (cio + /3o)ci2 
Therefore, we arrive at the photon field renormahzation constant 



ZA(a) 



4 a 



-2e 



a \' 



+ 



3 4:718 V47r£- 
and the renormahzed photon propagator at fi^ = —p^ (and e = 0) 



-P 



9 Att 



+ I6C3 



55\ / « \2 



The anomalous dimension ()3.45j) of the photon field is 



7a 



dlogZ 



A 



d 



d log /i d log n 

a 



a 
4:7ie 



+ ^20 - 7:^1 + Z2ie 



a Y 
Aire J 



4ne 



(-2^-2/5o£)+2(.2o-^^? + 



a 



2^20 ~ 772^1 + t:PoZi ) — 1-221 



a 
Aire 



(5.22) 



(5.23) 



(5.24) 



(5.25) 



(5.26) 



(5.27) 



It must be finite at £ ^ 0. Therefore, Z20, the coefficient of in the two-loop term in 
Za, cannot be arbitrary. It must satisfy 



Z20 



(5.28) 



In other words, C20, the coefficient of in the two- loop term in DxijP') ()5.19|) . must 
satisfy 

C20 = 7tCio(cio + l3o) ■ (5.29) 



Then 



lAO 



-Az. 
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i.e., the coefficients in the anomalous dimension are determined by the coefficients of 1/e 
in Za- Therefore, Za_ must have the form 



1 a 1 / \ 2 

Za = 1- ipAo^^ + g [lAoilAO + 2/5o) - 27Ai£] j + ■ ■ 



(5.30) 
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This is indeed so (see ()3.89|) for /?o), and we obtain 



This can also be understood in a shghtly different way. The information contained 
in Za is equivalent to that in 7^(a). This renormalization constant is gauge- invariant, 
because n(p^) is gauge-invariant in QED (there are no off-shell external charged particles 
in it; this is not so in QCD, where gluons are "charged"). Therefore, Za depends on /i 
only via a(yu) (there is no a(/i) in it). Dividing ()3.45j) by ()3.85j) . we obtain 

rflog^A _ 1 7A(a) _ 7A(a) /?(Q)7A(a) 

dloga 2e + /5(a) 2e "' ^ ' ' 

Any minimal ()3.8|) renormalization constant can be represented as 

ZA = exp(^^ + ^ + ---^ , (5.33) 

where Zi starts from the order a, Z2 — from a^, and so on. Then 

dZi 1 dZ2 1 ^/ \ / \ 

o7A(a) , — = -/?(a)7A(a) , . . . 



dloga 2 dloga 2 

and 

^ 1 f"' , .da 1 a 1 /a\^ 

1 r da 1 ^534) 



One can obtain 7a(«) from Zi, the coefficient of 1/e in Za, and vice versa. Higher poles 
(1/e^, . . . ) contain no new information: at each order in a, their coefficients {Z2, Z3, 

. . . ) can be reconstructed from lower-loop results. Up to two loops, this gives us ()5.3()j) . 

Now let's return to ()5.20|) with arbitrary L. It should be equal to the product of ()5.21|) 
and ()5.22|) . Equating a terms, we obtain 

Zi = Cio, ri{L) = cn - CloL , ru{L) = Cu - CuL + Ciq— ■ 

Equating terms, we obtain 

Z20 = C20 - PoCio , Z21 = C21 - 2c2oL - (cio + /3o)(cii - ciqL) , 
r2{L) = C22 - 2c2iL + 2c2oL^ - (cio + Po) ( C12 - cnL + cio^ ) . 
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But the renormalization constant Za cannot depend on kinematics of a specific process, 
i.e., on L. Therefore, terms with L in 221 niust canceL This is ensured by the consistency 
condition ()5.29j) . The renormahzation constant is given by ()5.30j) . and the renormahzed 
propagator at £ = — by 



a 

+ '4^ 



r2(0) + i (7A1 + ri(0) (7^0 - 2/3o)) L + ^7ao(7ao - 2/?o)i^' 



(5.35) 



+ 



This resuh can be also obtained by solving the RG equation ()3.47p . Taking into account 

we find that the coefficients obey the following equations 

drAL) 1 dr^iL) 1 (\ ^\ 
o7o, ^V^=o7i+U7o-/?o 



dh 2' ' dL 2' V2 



Solving them, we reproduce ()5.35p . 

Substituting the coefficients, we obtain 



.=.,,V)^l.t^(i.-f)H-(t^) (^-i..lOC,-f).... ,5.3, 

This is a typical example of perturbative series for a quantity with a single energy scale 
(— p^). In principle, we can choose the renormalization scale [i arbitrarily: physical results 
don't depend on it. However, if we choose it far away from the energy scale of the process, 
the coefficients in the series contain powers of the large logarithm L, and hence are large. 
Therefore, truncating the series after some term produces large errors. It is better to choose 
the renormalization scale of order of the characteristic energy scale, then \L\ < 1. The 
coefficients contain no large logarithm and are just numbers (one hopes, of order 1). The 
convergence is better. 

For example, when describing QED processes at LEP, at energies ~ my^, it would be 
a very poor idea to use the low-energy a with ~ mg: coefficients of perturbative series 
would contain powers of a huge logarithm log(mM//me). Using a{mw)-, which is about 7% 
larger, makes the behaviour of perturbative series much better. 

This renormalization-group improvement works for all processes with a single character- 
istic energy scale. If there are several widely separated scales, no universal method exists. 
For some classes of processes, there are some specific methods, such as, e.g., factorization. 
But this depends on the process under consideration very much. 
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5.4 Charge renormalization 

In QED, it is enough to know Za to obtain charge renormahzation ()3.77|1 . With two-loop 
accuracy, we have from ()5.25l 

1 4a/ 16 , ^ 



The /5-function ()3.86|) is simply 



2 

or (see jOH)) 



/3(a) = ~\lA{a) , (5.38) 



/?(«) = -1^- 4 f^) +■■■ (5.39) 



3 47r V47r 

The information contained in is equivalent to that in /5(a). Dividing ()3.86|1 by ()3.8f 
we have 

^log^g /3(a) /3(a) /52(a) 

ct log a e + /5(a) e e"' 
Writing Z^ in the form ()5.33p . we have 

dZx dZ2 \ 

-/5(a), -; = 13 (a),... 



d log a (i log a 

and 



a y 



'0 

° 2,„A^ _ 1^2 / iiV^ (5.41) 







One can obtain (3{q) from Zi, the coefficient of 1/e in Zq,, and vice versa. Higher poles 
(l/e^, 1/e^, . . . ) contain no new information: at each order in a, their coefficients (Z2, Z3, 
. . . ) can be reconstructed from lower-loop results. Up to two loops, 

^°=^-*£+('^«-H (£)'+■■■ ("2' 

Comparing this with ()5.37j) . we again obtain ()5.39j) . 

We shall later need /?(a) for QED with n/ massless lepton fields (e.g., electrons and 
muons at energies ^ m^). The photon self-energy at one loop ()3.34|) fFig. irUj) and at two 
loops ()5.18p (Fig. UH) gets the factor n/ (at three loops, there are both n/ and terms). 
If 1 — n(p2) is expressed via the renormalized a(/i), it should be equal to 

l-Il{p') = Z^{a{^^))[p'Dl{p'■^i 
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where Zc,{a) is minimal ()3.8p . and the second factor is finite at e — 0. Now we shall follow 
the derivation in Sect. 15.31 At one loop ()5.23p . zi gets the factor nj, because cio gets it 
(from ni(p^)). Therefore, /5o becomes nj times larger. At two loops ()5.24j) . Z21 also simply 
gets the factor Uf from C21, because cio + Po = zi + Po = for ()5.42j) . Therefore, 

/3(") = -|«/:^-4",(^)%-- (5.43) 

At three loops, both n/ and nj terms appear. The photon field anomalous dimension is, 
from (jSSHl), 

= + (-) +■•• (5.44) 

It is easy to write down the renormalization constants Za fj5.42|) and Za ()5.30|) . 

5.5 Electron self-energy 

The electron self-energy at two loops is given by 3 diagrams (Fig. 



Figure 44: Two-loop electron self-energy 




k + p 



Figure 45: Photon self-energy insertion 



The first diagram (Fig. contains one-loop photon self-energy subdiagram 

i {k'^g^iu — k^ky) n(A;^). This tensor is transverse, and hence longitudinal parts of the pho- 
ton propagators drop out. Therefore, this diagram is gauge-invariant: 

-illvii) = j (^^eo7^^ (f^^p ^eo7'^ * i^'^g^, - k^k,) . 

We shall set = — 1; the power of —p^ will be restored in the result by dimensionality. 



54 



Taking ^ Tr|i of both sides and substituting the result ()3.34p . we obtain 



el d-2„ [ d'^k N 



(47r)^/2ci-l 'J {2'kYDiD 
iV = 1 Tr |^7'^(^ + ■ {k'g,. - KK) . 

Using the "muhiphcation table" (j3.32p and omitting Di in the numerator, we have 

N^{d- 2)Dl -{d~ 3)D2 - 1 , 

and ^ 

^v-i = TT^T^Gi [Gil, 2 + e) + id- 3)G(1, l + 6)~{d- 2)G(1, e)] . 

Finally, using ()3.59p and restoring the power of — p^, we arrive at 




Figure 46: Electron self-energy insertion 

The second diagram (Fig. contains one-loop electron self-energy subdiagram 
— iSy((A; + p)'^){^ + 'i>)- For simplicity, we shall calculate it in the Feynman gauge oq = 1: 



Taking | Tr|^ (with = —1) and substituting the result (|3.6ip . we obtain 

ei d-2^ [ d'^k N 1 ^ , 

Using (j3.32j) and omitting D2 in the numerator, we have N ^ {d — 2){Di + 1), or 

Ev^2 = J^,^-^-^G^ [G(l, l + e) + G{l, e)] . 
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Finally, we arrive at 

- • ^^-^^^ 

The result in the arbitrary covariant gauge is very simple: 

eg(-p-)-'g^(^-2»^-3)^.^ (5.47) 



J 1/2 



{Any d-A 



The reason is following. It is not difficult to show (using ()3.59p ) that the diagram of 
Fig. ^1 (see ()3.56|) ) with the denominator Di raised to an arbitrary power n instead of 1 is 
proportional to Oq (prove this!). This means that in order to calculate this diagram with 
an arbitrary insertion(s) into the electron line, we may take the upper photon propagator 
in the Feynman gauge, and then multiply the result by oq. 




Figure 47: Truly two-loop diagram 



The third diagram (Fig. ITfj) is truly two-loop; it has the topology of Fig. |221 In the 

Feynman gauge. 



d'^ki d'^k2 
(2^(2^ 



(ki - k2y''''^^'\k,+py^^'' ki {k2 + pf 

d'^ki d'^k2 N 
{2nY{2nY D,D2DsD^D,' 



or 



4 



N = -TT^^^2Y{h - hh,ih+^hu. 

Using the "multiplication table" ()5.15|) and omitting all products DiDj except -D1-D4 and 
we have 



d-2 



{d - A)DiDi -{d- 8)D2D3 



+ 2(Di -D2-D3 + Di) + 2{Dl + Dl) + {d- 4) A 
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Therefore, 



4G(0, 1, 1, 0, 1) + 4G(1, -1, 1, 1, 1) + - 4)G(1, 1, 1, 1, 0) 



What is ^(1, —1, 1, 1, 1)? It would be trivial to calculate first the inner loop integral 
{cfik2) and then the outer one {cfiki), if there were no D2 in the numerator. In principle, it 
is possible to calculate such integrals using integration-by-parts recurrence relations. But 
we have not discussed the necessary methods. Therefore, we shall calculate this integral 
in a straightforward manner: 

^, , If d'^ki f d% r ,21 



1 r d^k, r d^k2 . , . 



d% 



D^D5 

In the second term, the inner vector integral is directed along ki, and we may substitute 

k2-k, 1 
k2 —^'^^ ^ 2 

(the terms with D4 or in the numerator give vanishing integrals). Therefore, 



J D,Ds2^ ''J 



d'^k2 



[G'(1,0,1,1,1)-G'(1,0,0,1,1)] . 



Using the trivial integrals 



Q 

Gil, 0, 0, 1, 1) = 6-2 , ^(1, 0, 1, 1, 1) = -^^G2 



d-A 



we obtain 



Substituting the integrals, we obtain the Feynman-gauge result 



e^i-p^)-^' d - 2 



(47r) 



(d-A)Gl + 8^G, 



d-A 



(5.48) 



(5.49) 



A similar, but more lengthy calculation in the arbitrary covariant gauge yields 



eti-pT''d-2 
(AttY a 



{{d - 2)al + d- 6)Gl - ^^r-^dd - A)al -d-A) 

(Jj "rr 



(5.50) 
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We shall consider a more general case of QED with n f massless lepton fields; the results 
for the usual QED can be easily obtained by setting nf = 1. Then the diagram of Fig. 1^ 
gives ()5.45|1 with the extra factor n^; the diagram of Fig. IHl gives ()5.47|1 . and that of Fig.lTfl 
gives ()5.5()j) . The complete result for the two-loop term in Sv'(p^) is 



-P 



2\-2e 



(An) 



-{d-2) 



2'L^G2nf -\{{d- 2)al + - 6) 



d-6 
ld-3 

+ 



2d- A 



{{M - S)al - d - A)G2 



(5.51) 



5.6 Electron field renormalization 

Now we proceed as in Sect. 15.31 The electron propagator ^S{p) has the form similar 
to (j5T9|l . but now the coefficients are gauge- dependent: 

From (jTFT]l and using fTI^i and we have 



e'f2{e) 



1 - 2e 



2e{l-e) 



_(l + £)(l-2e)(l-3£)(2-3e) 



g2nf 



1 + e 2 



[1 - 2ey 



4 - £ - (2 - ?,e)al 
2(l-3e)(2-3e) 
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We can expand these functions in e: 



e^e^'^^/2(£) = C20 + C2ie + C22e^ H 

In our particular case, /{(e) = 0, 

1 



-l-e 



e'e^'^'f2{e) = | + (n, + + ^) . + Qn, + 3ag - \C2al + ^ 



The situation here is simpler than in Sect. 15.31 the terms with gf and with g2 in e'^f2{e) 
are separately finite at e — 0, so that we can put g2 = gf with the e"^ accuracy, and Cs 
does not appear. 

The propagator expressed via the renormalized quantities «(//), a(/x) is, with the 
accuracy, 

^Sip) = 1 + ^e-^^Z„5e^^ [f[{e) + f^is)ZAaifi)] + [£-^\-^^^e'e'^^ f2{e) , 
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see dSIini)- At L = we have 



+ 



a \' 



4ne 



) 



C20 + C2ie + C22e^ H Po (cio + cne + c^e^ H ) 



+ ^lAoa (c'/o + c'lV + c'l^e^ + ■■■) 



Here ci„ = c']^„ + c'l^a^fi); in C2„, we may substitute a(/i) instead of qq. This should be 
equal to Z^{a{fi),a{fi));^Sr{p; fJ^). Equating a terms, we obtain (j5.23|) : equating terms, 
we now have 

2^20 = C20 — PqCio — -'jAoCioa , 

2:21 = C21 - (cio + Po)cu - ^jAoc'lia , 

r2 = C22 - (cio + /5o)ci2 - ^7AoCi2a • 

There are extra terms as compared to the gauge- invariant case ()5.24j) . In QED 7^0 = —'^Po- 
We obtain 



Z^{a, a) = 1 — a- 



a 
'Aire 



+ 



4 



Arce 



) 



+ 



(5.52) 



It is also easy to write down fSr{j>] /i) at L = (i.e., /x^ = — p^), or even at arbitrary /i, 
but we shall not do this here. 

The anomalous dimension of the electron field is 



7v 



d log 



a 

-2zi 4 

477 



d d 
-2(£ + /3(a))Trr— -7A(a^ 



9 log a 



(9 log a 



[z[ + z'la) 



a 



1 



a 



Z20 - 0% + 0/^0^1 + -TlAoZia - + 2:21 



4 



a y 



where zi = z[ + z'la. It must be finite at e — >^ 0; therefore, the two- loop term must satisfy 
the self-consistency relation 



or 



2^20 = ^^1(^1 - /^o) - ^lAQz'la . 



1 1 /, 

C20 = -^cioicio + (3o) + --fAoc'loa . 



(5.53) 
(5.54) 
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The renormahzation constant of a non-gauge-invariant quantity (such as the electron field) 
must have the form, up to two loops. 



1 Q; 1 / \ ^ 

Z^ = l- 2^^o^ + g [li,o{.l^o + Wo) + lAolloa - 27^15] (^^j + ■ • • (5.55) 

where 7^0 = 7^o + I'iio'^ (this formula generalizes the gauge- invariant case ()5.30p ). 
We see that the self-consistency condition is indeed satisfied, and 

7^(a, a) = 2a— - (4n; + 3) (— j + • ■ ■ (5.56) 

In the normal QED, with just one charged lepton field, 

7^(«,a) = 2a— -7(^— j +■■■ (5.57) 

The electron mass anomalous dimension can be found in a similar way. We can calculate 
T^sip^) at two loops (neglecting m^) by retaining mo in the numerator of a single electron 
propagator in Fig. and setting mo ^ in all the other places. This single mo has 
to be somewhere along the electron line which goes through all the diagrams, not in the 
electron loop in the first diagram: we need one helicity flip of the external electron, and 
one helicity flip in a loop yields zero contribution. Then we extract Zm from ()3.99p . It 
must be gauge-invariant, because m(/i) is gauge- invariant. The result is 



a/ 20 \/a\2 



In the normal QED, with just one charged lepton field. 



, , a 11 / a \2 , ^ 

= «S - T (4^) +■■■ f^'^^' 

5.7 Two-loop corrections in QCD 

The one-loop /3-function in QCD is ()4.45|) : the two-loop one is 

34 20 

A = —C\ - ACpTprif - —CATprif. (5.60) 

The second term follows from the QED result ()5.39p : non-abelian terms (containing Ca) 
are more difficult to derive. 

How to express a<((/i') via as(/i), if the scales /i and fi' are not too widely separated? 
We need to solve the RG equation ()4.46j) with the initial condition at /i. Let's introduce 
short notation: 



as(/i) , _ a^(/i') , _ A 
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Then the integral of our equation is 



1 + bias H al pi 

It is natural to introduce notation 

l = 2/3o\og^. 
Expanding the left-hand side, we obtain 



' '"'-20,10^^. (5,61) 



1 - ha, + ...)^ = l- l + 6,log^^ + ... = i. (5.62) 



The solution is a series: 



a[ = a.. 



(l + Citts + C2a1 H ) . 

Substituting it into ()5.62|) . we obtain 

1 1 

— [l - citts + {cl - C2)a1 H 1 \-bi [cias H ] = / 

as as 

Therefore, ci = —I and C2 = — bil. The final result is 

[1 - las + {l^ - bil)al + ■ ■ ■] . (5.63) 



a'g = as 



It is easy to find a^, ... corrections; they will contain 62, • • • 

It is not possible to find the exact solution of the RG equation for as{fi) in elementary 
functions, if two (or more) terms in (3{as) are kept^. However, an implicit solution can be 
obtained. Separating variables in the RG equation ()4.46|1 . we can write it as 

= — dlog/i . (5.64) 



2(3{as) as 

Let's subtract and add two first terms in the expansion of the integrand in a^: 

Wias) 2Poas^2Pi)as 2/5oa.(/i) 2/?^^ " A^s ' ^ ^ 

The added terms are explicitly integrated; the difference behaves well at Oa — 0, and can 
be integrated from 0. The integration constant has appeared here, as in the one-loop 
case fl4.47j) . We can solve for this constant: 



A = /iexp ( ) a,(/x)-'^^/(2'^o)ir(a,(/i)) , (5.66) 



^With two terms, the solution can be written via the Lambert ly-function. 
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where 

1 f3i \ da 



/f(a,)=expy^ [^—- — + ^J- = 1 + ... (5,67) 

This is a regular series in a^; its term contains (32, and so on. 
The quark mass anomalous dimension is, up to two loops, 



7r> 



6C,^ + C, (3C. + f - f T,n,) (^)V • • • (5.68) 



Everything here, except the Ca term in 7^1, can be obtained from QED result ()5.58|) by 
inserting obvious colour factors. The solution ()3.109|) of the RG equation can be written 
in the following convenient form. Let's subtract and add the first term of the expansion of 
the integrand in a^; the difference can be integrated from 0: 



7mo/(2/3o) 

Q;(/i) 



m(/) = m(/x) ( ^ ) K,^{a,{^,'))K-^{a,{^^)) , (5.69) 



where for any anomalous dimension 



2 

+ ■ 



we define 

7("s) 7o \ das 



K^{as) = exp 



Ols 



\2(3{as) 213 J a 



(5.70) 



_ ^ 7o / 71 f3i\as 
2/5o V70 /5o;47r 

This function has the obvious properties 

Koi^as) = 1 , K^j{as) = -ft'^^(as) , K^^+j^ios) = K^^{as)Ky,{as) ■ 
The solution ()5.69|1 can also be rewritten as 

m(;x) = m(^^j K,^{as{^^)), (5.71) 

where m is a renormalization group invariant which characterizes the given quark flavour. 
The running mass m(/i) decreases when /x increases, because 7^0 = ^Cp > 0, /5o > 0, so 
that the exponent in ()5.71|) is positive, and a<j(/x) decreases with /i. 

The QCD [3- and 7-functions are also known at three [131 CH HSj and four [121 El UHl 
[inil2Dlini loops. 
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6 On-shell renormalization scheme 



6.1 On-shell renormalization of photon field 

We shall consider QED with non-zero electron mass ()3.94|1 . Until now, we used MS renor- 
malization scheme. In this scheme, the electron mass m(/i) and the coupling depend 
on the renormalization scale /i. It would be more exact to call MS a one-parameter fam- 
ily of renormalization schemes. But what are the experimentally measured electron mass 
m ~ 0.511 MeV and coupling a ~ 1/137? They are defined in another renormalization 
scheme — on-shell scheme. It contains no parameters. MS scheme is most useful at high 
energies, when the electron mass can be neglected (or considered as a small correction); /i 
should be of order of characteristic energies. When energies are of order m, the on-shell 
scheme is often more convenient. We shall discuss it now. 

The photon field renormalized in the on-shell scheme A^^ is related to the bare one Aq 

by 

Ao = {Z°^f'A,,, (6.1) 

where Z*^ is the renormalization constant (it is not minimal in the sense fl3.8p ). The 
renormalized propagator is related to the bare one by 

D^ip') = Z°^D°^{p') . (6.2) 

The bare photon propagator near the mass shell is 

We require, by definition, that it behaves as the free propagator 1/p^ near the mass shell. 
Therefore, the photon-field renormalization constant in the on-shell scheme is 

ZT = . (6.4) 

^ i-n(o) ^ ' 

How can we calculate 11(0)? From ()3.15|) we have 

U^.ip) = {d - l)p'U{p') , 

and hence 

2c/(c/ - l)n(O) . (6.5) 



dpu dp 

The left-hand side has zero external momenta — it is a vacuum diagram, with the only 
energy scale m. It is not difficult to calculate a few more terms of expansion of n(p^) in 
p^ in the same manner. 

At one loop (Fig. I48|l. setting tjiq = 1 (it will be restored by dimensionality), we have 



d d 
dpy dp^ 



So{k + p) 



(6.6) 

p=0 
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k + p 



Figure 48: One- loop photon self-energy 



Expanding the free electron propagator 5*0 (A; + p) (j3.95|) in p, 



So{k+p) 



+ 1 



+ 



^ 2p-k{^ + r 



fc2-l k^-1 (A;2-l)2 

2p-k^ p\^ + l) {2p-kf{^+l) 



we obtain 



d d 
dpu dp'' 



{k^ - 1)2 (P - 1)2 



+ 



Soik + p) 



- 1)3 

8k\^ + l) 



+ 0{p') 



p=0 



(P - 1)2 



+ 



Using also (j3.27|l . (j3.2fij) to simplify 'y^So{k)jf^, we can easily calculate the trace: 



n(o) 



d{d -I) J (27r) 



d'^k 



4(c/-3) d^-Ad + A 



D2 



-4 



(4vr)'^/2 



[8^(4) + A{d - 3)V(3) + {d^ -Ad + A)V{2)\ 



where D = 1 — k'^, and the definition ()2.1|) of the vacuum integrals was used. Using ()2.7p . 
we can express all V{n) via V{2) = T{e); also restoring the power of niQ by dimensionality, 
we finally arrive at 

m) = -f^ne). (6.7) 

The ultraviolet divergence of the diagram of Fig. [TUldoes not depend on masses and external 
momenta, and is the same in (j6.7p as in (j3.35j) . 

Therefore, the photon-field renormalization constant in the on-shell scheme is, with 
one-loop accuracy. 



'70s 



1 - 



A elm, 



"0"''0 

3 {ATiy/"^ 



Tie) 



(6. 



The renormalized photon propagator in the MS scheme and in the on-shell scheme are 
both ultraviolet-finite, therefore, the ratio Z'^ / Z must be finite at e — > if both 
of them are expressed via renormalized quantities (the on-shell photon propagator has 
no IR divergences because the photon is neutral). Therefore, we can reproduce the MS 
renormalization constant ()3.42j) from ()6.8|) . 



64 



6.2 One-loop massive on-shell propagator diagram 

Before discussing the electron mass and field renormalization in the on-shell scheme, we 
have to learn how to calculate the relevant diagrams. 




Figure 49: One- loop on-shell propagator diagram 
Let's consider the on-shell propagator integral (Fig. HHl = "^^) 



ni,n2) 



D^=m^-{k + pf = -k'^-2p-k, D2 = ~P 



(6.9) 



The power of m is evident from the dimensional counting, and our aim is to find the 
dimensionless function M(?t,i, -0.2); we can put m = 1 to simplify the calculation. It vanishes 
if rii is a non-positive integer. 

Using Wick rotation and a parametrization (j2.3j) . we rewrite the definition (|6.9|) as 



M{ni,n2) 



'd/2 



(6.10) 



We want to separate a full square in the exponent; to this end, we shift the integration 
momentum: 



k' = k + 



ai + 0L2 



p 



and obtain 
M(ni,n2) 



TT 



-d/2 



r(ni)r(n2) 
1 







1 exp 


ai + a2_ 










al 


1 exp 


ai + a2_ 





T{ni)T{n2) 

Substituting t] = ai + a2, «! = rjx, a2 = ri{l — x), we get 



(ai + a2)-'^^^a'l'-^a^^-^dai da2 . 



M(n.,.2)^^zr^/^--^(l-r-^^-f 



Tin,)] 

T{-d/2 + ni + n2) 



T{ni)T{n2) Jo 



X 



d—ni—2n2 — l i 



g-'?a;^^-d/2+?ii+n2-l^^ 



1 -x)"2-Mx. 
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Finally, we arrive at 



, . r -rf/2 + ni + ^2 r - - 2^2 , . 

M m, ^2 = ^ . 6.11 

r(ni)r(d - ni - n2) 

The denominator in (j6.9p behaves as (/c2-jni+n2 at A; — oo. Therefore, the integral 
diverges ii d> 2(ni + ^2). At ci — > 4 this means ni + n2 < 2. This ultraviolet divergence 
shows itself as a 1/e pole of the first F function in the numerator of 1)6.111) (this F function 
depends on rii + n2, i.e., on the behaviour of the integrand at /c — 00). The integral ()6.9|1 
can also have infrared divergences. Its denominator behaves as /c"i+2"2 _j. g^^^j 
integral diverges in this region ii d < rii + 2n2. At d ^ 4 this means rii + 2n2 > 4. This 
infrared divergence shows itself as a 1/e pole of the second F function in the numerator 
of ()6.11|) (this F function depends on rii + 2n2, i.e., on the behaviour of the integrand at 



Figure 50: One- loop on-shell propagator diagram 



Q 

Figure 51: The basis integral 



Let's summarize. There is one generic topology of one-loop massive on-shell propagator 
diagrams in QED and QCD (Fig. IHU)) . All Feynman integrals of this class, with any 
integer indices ni, n2, are proportional to Vi ()2.8)1 (Fig. lHT)) . with coefficients being rational 
functions of d. For example, for M(l, 1), 



Two-loop [231211122] and three- loop |2SII2E] massive on-shell diagrams can be calculated 
using integration by parts. 



6.3 On-shell renormalization of electron mass and field 

The electron mass m in the on-shell renormalization scheme is defined as the position of 
the pole of the electron propagator. On-shell external electron lines have = m^; it is 
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convenient to use the free propagator containing the same quantity m rather than itlq. 
Therefore, we rewrite the Lagrangian ()3.94|) as 

L = Tpoiilp — m)ipo + STmpQipo , 6m = m — mo , (6.13) 

and consider the mass counterterm not as a part of the unperturbed Lagrangian, but as a 
perturbation. Then the free electron propagator 

Soip) = = 4^ (6-14) 
p — m — 

contains m, and the mass counterterm produces the vertex 

— ^ — x-^ — =i5m. (6.15) 
The mass renormahzation constant in the on-shell scheme is defined as usuah 

mo = Z^m. (6.16) 
It is more convenient to write the electron self-energy ()3.97j) in the form 

S(p) = mJ:i{p^) + m)S2(p') (6.17) 
now. The electron propagator 

^^^^ " [1 - ^2{j>^W -m) + 5m- mSi(p2) ^^'^^^ 
has a pole at = m^ if 

(5m = mSi(m^), (6.19) 

or 

OTfl 

ZZ = \ = 1 - Si(m2) . (6.20) 

m 

The equation ()6.19j] can be solved for bm by iterations (at higher orders, its right-hand 
side contains bm because of the vertex ()6.15j) . 
Near the mass shell, we can expand Si(j9^) as 

Si(p') - — = S;(m2)(p2 _ + . . . (6.21) 
m 

so that 

1 — h2[m'^) — 2m^L^(m"^) p"^ — m^ 
The electron field renormalized in the on-shell scheme, 

^0 = (^;^) i^os , S{p) = Z;^SUp) , (6.23) 
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is defined in such a way that its propagator Sos{p) behaves as So{p) ()6.14|) near the mass 
sheU: 



= [1 - S2(m=^) - 2m^J:[{m^)] ' . (6.24) 



In order to calculate and Z^*^, it is convenient to introduce the function 
T{t) = ^ Tr(^ + l)S(mt;(l + t)) = ^i{m'^) + [Salm^) + 2m'^J:[{m'^)] t + ■ ■ ■ (6.25) 



so that 



z- = i-T(o), z;^ = [i-T'io)r' 



(6.26) 




Figure 52: One-loop electron self-energy 

Let's calculate it at one loop fFig. I52|). We put m = 1; the power of m will be restored 
by dimensionality: 

d'^k 1 1 



T(t) = -let 



{2nY D,{t)D2 



i Tr(^ + 1)7'^(|^ + + l)r (^9,. + , (6.27) 



where 



p = vil + t), Di{t) = 1 - {k + p)\ D2 = -P. 
In calculating the numerator, we can express 

p-k = ^[D2-D,{t) + l-il+tf] 

and omit terms with Di{t), because the resulting integrals contain no scale. We obtain 



T(t) = -tel 



d'^k 1 



2 d- 2 



This result is gauge-independent. 



i-t) + o{e 



Expanding 



Di{t) =Di + {Di - D2 - 2)t + 0{t^) , Di = l-{k + v) 



we obtain 



T{t) = -lel 



d'^k 
{27iY 



2(1 -t) (rf- 2)(1 - 2t) _^ At {d-A)t {d-2)D2t 



D1D2 



2Di 



DID2 Dl 



2Dj 



68 



Using ()6.1H) and restoring the power of m, we finally arrive at 



i-W = ^r(.)^(i-i) + 0(t^). (6.28) 



zS = z? = i-Si7^r(£)5-4- (6-29) 



Therefore, 

elm-^' d-l 
(47r)'^/2 ^>d-'i 

The fact that = Z'!^ at one loop is a pure accident; at two loops it is no longer so [2S1I21!- 
The on-shell mass m is a measurable quantity, and hence gauge- invariant; therefore, Z'^ is 
gauge-invariant to all orders. It has been proved that Z'!^ is also gauge-invariant in QED 
(but not in QCD). 

What's the relation between the on-shell mass m and the MS mass m(/i)? 



mo = Z^{a{fi))m{fi) = Z^m , 

and therefore 

m{fi)= „ ... m. (6.30) 

We have to re-express Z!^ via a(/i); then the ratio of renormalization constants is finite at 
e ^ 0, because both renormalized masses are finite. At one loop 

"^^^-l-6f flogii + ^'l . (6.31) 



m Att \ m ^ 

This formula is OK when fi m; otherwise, it is much better to relate m{m) to m, and 
then to solve the RG equation with the initial condition m{m). 



6.4 On-shell charge 

On the mass shell (p^ = m^), r'^(p, p) has only one 7-matrix structure, 

r^(p,p) = Z{JV, (6.32) 

if sandwiched between U2- ■ - Ui which satisfy the Dirac equation (see Sect. 16.51 for more 
details). The physical matrix element of scattering of an electron by a photon is 

e.V^^Z'^iZ^lf' . 

The only case when all 3 particles are on-shell is = m?, and the photon momentum 
g — >■ 0. By definition, this matrix element is 

e7^ (6.33) 

where e is the renormalized charge in the on-shell scheme. It is related to the bare charge 
by 

e, = {Z':f'\. (6.34) 
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Therefore, 



(6.35) 



On-shell electron charge is measured in macroscopic experiments with smooth electromag- 
netic fields (having g — > 0); it is 

1 

a = — ^ . (6.36) 

ATI 137 ^ ' 

In QED, the situation is simplified by the Ward identity ()3.75|) . Near the mass shell 



S{p) 



yos 



and hence 
Therefore, 
and 

At one loop, from ()6.8j] . 



^ — m ' 



Z^'Z^' = 1 



1 + 



:r(e] 



3 (47r)'^/2 

The MS (running) coupling Q;(/i) is related to the on-shell coupling a by 

1 



(6.37) 
(6.38) 

(6.39) 



ZaW 



4a(^///.x- 
3 Att \\mJ ^ ' e 



Therefore 



a(/i) = a 



8 a , a 
1 + log^ 
SAtt m 



(6.40) 



We can always find the /x-dependence from the RG equation, the initial condition is a{m) = 
a at one loop. 

6.5 Magnetic moment 

Let's consider scattering of an on-shell electron in electromagnetic field. The physical 
scattering amplitude is 

eoZ^^Z^^^f' u'T'^ip,p')u, (6.41) 

where = m^, p'^ = m^, and the initial and final electron wave functions satisfy the Dirac 
equation 

-^u = mu , u'f = u'm . (6.42) 
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We can substitute the on-shell charge 

e = e,{Zff'\ 

Using the Ward identity ()3.74|) and the electron self-energy ()6.17p . we obtain 

where q = p' — p, and hence for on-shell electron wave functions ()6.42|) 

u'Vq^u = . (6.43) 

What 7-matrix structures can u'V^u have? Using the Dirac equations ()6.42|) . we can 
always eliminate |^ by anticommuting it to m, and f by anticommuting it to u'. For example, 

n^'a^'^q^u = i u' {{p + p'f - 2m-i^) u . (6.44) 

We are left with 3 structures: 7^^, {p + p'Y, and q^ = {p' — p)^. The last one is excluded 
by ()6.43|1 . Therefore, 



u = u' 



{FM') + F,{q^))l'-F2{q')^^^^'^' 



2m 



u . 



Using ()6.44p . we can also rewrite this as 



(6.45) 



u = u' 



2m 



u 



u' 



2m 



mq') + F,{q^))^ 



u 



Let's rewrite ()6.45|1 as 



{p + p'Y 

2m 



then 



Fi(g2) = Z-(A + /2), F,iq') = -Z^h. 
If we introduce the traces 



then 



4m^ 



—^TTT^f + m)T^^{^ + m) 
and we can find fi by solving the linear system: 



(6.46) 



(6.47) 
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Let's introduce the notation 



4m^ 



Calculating the traces, we find 

M = 2 

and hence 



[l + tf l + t 
l + t l-{d-2)t 



1 



l-{d-2)t -(l + t) 

2(c/-2)t(l + t)2 \^ -(l+t) (1 + t)^ 



Finally, we obtain from ()(i.47j) 

1 



^■0 



2(rf-2)(l +t)2 4 
'70s 



Tr [{d - l)v, - (1 + th,] if + l)r'^(^ + 1) 



(6.48) 



1 



2(rf - 2)t(l + ty 4 [(^ - - 2)t)t;, - (1 + th,] if + l)r'^(^ + 1) , (6.49) 

where v = p/m, v' = p' /m. At the tree level, = 7^^, and, naturally, we obtain Fi{q^) = 1, 
-^2(1?^) = 0. When calculating loop corrections, we can apply these projectors to integrands 
of vertex diagrams, and express Fi^2{q^) via scalar integrals. 
The Dirac form factor (|6.48|) at = is 

1 



Fi(0) = Z;^ - i:iv^T^'{mv,mv){i) + 1) . 



Due to the Ward identity (pTT^ . 

d 1 



Fi(0) 



dpi^A 



TrE(p)(^ + l) 



p=mv 



ZT [1 - r 



(6.50) 



(6.51) 



see fj6.25|) . ()6.26p . The total charge of electron is not changed by radiative corrections. 

It may seem that the Pauli form factor ()6.49|) is singular at — > 0. Of course, it is 
not. Let's substitute the expansion 



Virnv, mv + q) = T^ + Ff ^ + 



m 



into ()6.49p . split q = {q ■ v)v + g_L (where q ■ v/m = 2t, q\/ 



m 



-4t(l + t)), and average 



over the directions of q± in the {d — l)-dimensional subspace orthogonal to v: 

At 



q- 



m 



2tv' 



d 



- [(l+t)(7"^-(l + rft)t;V] 



We obtain 

^2(0) 



'70s 

d-2 



■TT{j,-dvM^ + l) 



2 1 



(6.52) 
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This means that in order to calculate the anomalous magnetic moment we need the vertex 
and its first derivative in g at g = 0. 




Figure 53: One- loop anomalous magnetic moment 

As already mentioned, the tree diagram (Fq = 7^^, F^'^ = 0) does not contribute to 
^2(0). The first contributing diagram is one-loop (Fig. l53|) . All charged external lines are 
on-shell, therefore, this vertex diagram is gauge-invariant. We shall use Feynman gauge to 
simplify the calculation, and put m = 1 (the power of m will be restored by dimensionahty) . 
Then the one-loop vertex is 

We need to expand it in q up to the linear term. In this diagram, only one propagator 
depends on q. Using 

5(p + g) = g(p) + ^-^- f + 2p-g + --- (6.53) 

for p = k + mw, we obtain (m = 1) 

.2 f d'^k 7,(1^ + ^ + 1)7^(1^ + ^ + 1)7° 

d'^k 7„ [DiY + 2(A; + vy{^ + ?^ + 1)] + + 1)7" 



-'"V (2^ DfD, 

where Di = 1 — {k + vY, D2 = —k"^ (see fl6.9p ). We may replace Z'lp — > 1 in ()6.52p . because 
corrections are beyond our accuracy; the two contributions to F2{0) are 

iel r d'^k iVn 



d-2j {2TrY DID2 ' 

No = \ Tr(7^ - dv^)^^{l/^ + ^ + 1)7^(1^ + f + 1)7"(^ + 1) 
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and 



2iel 



{d-l){d 



{2nYDfD2 ' 



X [D,Y + 2ik + vYilt + ^ + 1)] 7^(1^ + + 1)7°(^ + 1) • 
When calculating A^o, we may omit terms with Df: 

No =^ -d{d - 2)DiD2 -{d~l){d- A)Di + ^d{d - 2)Dl + d{d - 3)L>2 - 4(rf - 1) . 

When calculating A^^i, we may omit terms with Df: 



3{d - 2)DfD2 + {d^ - 8d^ + 17d - U)Df - 3{d - 2)DiD^ 
- {d^ - 8d^ + 21d - 26)DiD2 + 4(rf - ifDi + (rf - 2)Dl + 4(rf - 3)Dl 
Collecting all this together, we obtain 



16A 



F2{0) 



1 



(47r)^/2 d-2 



d{d - 1)M(1, 0) - (d - l){d - 4)M(1, 1) 



+ -d{d - 2)M(2, -1) + d{d - 3)M(2, 0) - 4(rf - 1)M(2, 1) 



+ 



d-\ 



3{d - 2)M(1, 0) + {d^ - 8d^ + lid - 14)M(1, 1) 

- 3{d - 2)M{2, -1) - {d^ - 8d^ + 2ld - 26)M(2, 0) 
+ A{d - lfM{2, l) + {d- 2)M(3, -2) 

+ 4(rf - 3)M(3, -1) - 16M(3, 0) 



Using and restoring the power of m, we finally arrive at the anomalous magnetic 

moment 



i^2(0) 

As expected, it is finite at e — > 0: 



elm-^^^^^^ {d-A){d-h) ^ 



(47r)d/2 



^2(0) 



a 
2^ 



d-3 



+ 



(6.54) 



(6.55) 



As we can see from the second line in fl6.46|l . the total magnetic moment of electron is 

a 

Fi(0) + F2(0) = ! + — + •• • 

(in units of Bohr magneton). It is not difficult to calculate the two-loop correction to the 
magnetic moment, if one has a program implementing the integration- by-parts evaluation 
of two-loop on-shell propagator integrals. The three-loop calculation |23 was a major 
breakthrough. 
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6.6 Two-loop massive vacuum diagram 

Let's consider the two-loop massive vacuum diagram (Fig. l5l|l : 



The power of m is evident from the dimensional counting, and our aim is to find the 
dimensionless function V^(ni, n2, na); we can put m = 1 to simplify the calculation. It is 
symmetric with respect to rii ^ it vanishes if rii or n2 is a non-positive integer. Using 
Wick rotation and a parametrization ()2.3|) . we rewrite the definition ()6.56p as 

Vr / _ai(fc2+l)-a2(fci+l)-a3(fci-fc2)2 



= r(„,)rMr(n3)y^" (6.57) 



X a"' ^02^ ^a"^^ ^daida2da3d'^kid'^k2 . 




Figure 54: Two-loop massive vacuum diagram 



The (i- dimensional integral of the exponent of the quadratic form is the product of d 
one-dimensional integrals: 



H d 



[det A] 



d/2 ■ 



(6.5J 



This is the definition of such d- dimensional integration; note that the result contains d as 
a symbol. In our case, 

A = ( "1 + "3 -"3 
Y -as a2 + as 

and we obtain 

V{ni,n2,n3) 



T{ni)T{n2)T{ns) 
X {aia2 + ttiOs -|- a2a3)^'^^'^a^^~^a2^~^a^'^~^dai da2 da^ . 
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To calculate this integral, it is most convenient to choose the "radial" variable t] 
«! + ^2 (Sect. EH}, and substitute = rjXi. The integral in i] is trivial: 



X J {xiX2 + X1X3 + X2X3)~'^^'^ Xi^~^ X2'^~^ x^^~^ 5{xi + X2 — l)dxi dx2 dx3 

" r(ni)rK)r(n3) J ^ ^ij+^sj x, [L x,) X, dx.dxs 

Substituting X3 = xi(l — xi)y, we get 

T{-d + ni + n2 + n3) y^'^'^dy 



V{ni,n2,n-i) 



T{n,)T{n2)T{n,) J, {y + lyi^ 



1 

X / X 





The integral 

y^^-My V{n^)V{d/2-n^) 



/o {y + iyi^ T{d/2) 

is easily calculated using the substitution z = 1/ {y + 1)] the second integral is the Euler 
i?-function. We arrive at the result [7] 

. r(f -n3)r(ni +723- f)r(n2 + n3 -f)r(ni + 722 + ^3 -rf) 

V^(ni, n2, ^3) = T(d\rt — —7, • (6-59) 

r (|j V[ni)V[n2)T[ni + ^2 + 2n^ - d) 

This is the only class of two-loop diagrams for which a general formula for arbitrary Ui 
(not necessarily integer) is known. 




Figure 55: Two-loop vacuum diagram 

Let's summarize. There is one generic topology of two-loop vacuum diagrams in QED 
and QCD (Fig. 13^ . All Feynman integrals of this class, with any integer indices n^, are 
proportional to ()2.8|) fFig. with coefficients being rational functions of d. 

Three- loop massive vacuum diagrams can be calculated using integration by parts [22] ■ 
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Figure 56: The basis integral 



6.7 On-shell renormalization of photon field and charge at two 
loops 

Using the results of the previous Section, it is not (very) difficult to calculate the photon 
self-energy 11(0) at two loops. We apply ()6.5p to the diagrams of Fig. HH calculate the 
derivatives in p, and reduce the problem to the vacuum integrals ()6.56p . The result is 

2 eX"'^ 2,.^ (^-4)(5^^- 33^ + 34) 
= -3^4^^ did^) • ^^-^^^ 

Therefore, for the on-shell charge e we obtain 



-,2 



e 

-2e /I n I 7^ in^2 /^2^-2£ \2 



3 = = {z^^y = 1 - n(o) 



^ 3 (47r)'^/2 3^ (2 - e) (1 + 26) V (47r)^/2 > ^ 



(6.61) 



where n(0) at one ()6.7j) ) and two ()6.60|1 was used. In the one-loop term, we have to 
substitute eg = Z^^e^ and mo = Z^^m, with one-loop flO^ and Z^ (K^ . This 
results in 

^2 

f^°')"^ 



^0 



2 



4 2 -2£ l + 7^_4^3 / 2 -2e 

= 1--- — — Fe -4e- - — -tttF +■■■ 6.62 

3(471)^^/2 v; (2-e)(l -2e)(l + 2£) V(4vr)°'/2 ^ ^ ^ 

= ^ - 3 (S^F^^f^' - .(2 + 1& + . . . ) (^r(.) j + . . . 

On the other hand, for the MS charge we have 

(let's pretend for a moment that we don't know Za yet). The ratio e^(/i)/e^ must be finite 
at e 0. At one loop, this requirement gives zi = —4/3. Setting /i = m, we can substitute 



rF £ 



47re (47r)"'/2 Ane 
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in the two- loop term, because the differences are of a higher order (here a = e^/(47r)). 
Thus we obtain Z20 = 0, Z21 = —2. We have reproduced in MS fj5.37|) . ()5.25|) from 
our on-shell calculation. Finally, we arrive at the relation between the MS a{m) and the 
on-shell a: 



aim) = a 



(6.64) 



The three- loop correction has been calculated in [22j. If we need a{fi) for /i 7^ m, we can 
solve the RG equation ()3.90|) with this initial condition. 

6.8 On-shell renormalization in QCD 

QCD perturbation theory is only applicable at large momenta (or small distances). There- 
fore, it makes no sense to renormalize the gluon field, the coupling, as well as light-quark 
masses and fields, in the on-shell scheme. However, it is possible (and often convenient) to 
use this scheme for renormalizing heavy-quark masses and fields, at the same time using 
MS for a<((/i), the gluon field and light-quark masses and fields. One- loop renormalization 
constants can be trivially obtained from the QED results (j6.29j) : 

^- = ^Q=1-C^.^r(.)^. (6.65) 

Similarly, the relation between the MS mass m{m) and the on-shell mass m is, from ()6.31|] . 

!!^ = l-AC,^ + ... (6.66) 
m An 

The two-loop correction has been found in [22] and the three- loop one — in j^H]^- The 
quark- field renormalization constant Zq at two loops [21] is gauge-invariant, but at three 
loops j2n] not (unlike in QED). 



7 Decoupling of heavy-particle loops 
7.1 Photon field 

Let's consider QED with massless electrons and muons having mass M. When we consider 
processes with characteristic energies E <C M, the existence of muons is not important. 
Everything can be described by an effective low-energy theory, in which there are no muons. 
In other words, muons only exist in loops of size ~ 1/M; if we are interested in processes 
having characteristic distances much larger than 1/M, such loops can be replaced by local 
interactions of electrons and photons. 

The effective low-energy theory contains the light fields — electrons and photons. The 
Lagrangian of this theory, describing interactions of these fields at low energies, contains 

^It had been found numerically [27] before the analytical result was obtained. 
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all operators constructed from these fields which are allowed by the symmetries. Operators 
with dimensionalities di > 4 are multiplied by coefficients having negative dimensionali- 
ties; these coefficients contain l/M*^'"^. Therefore, this Lagrangian can be viewed as an 
expansion in 1/M: 

L' = ^P',^1P% - \F^,^Fr - ^ {d,KY + ^foF^^.^cr'^-'i^'o + O (7.1) 

(primed quantities are those in the effective theory). The coefficients in this Lagrangian 
(c, . . . ) are fixed by matching — equating S'-matrix elements up to some power of Pi/M. 

Let's imagine for a moment that we only know electrons and photons. There are two 
ways to search for "new physics" : 

• To raise energies of our accelerators in the hope to produce real new particles (e.g., 
muons) ; 

• To measure low-energy quantities (such as the electron magnetic moment) with a 
high precision in the hope to find effects of higher terms in the effective Lagrangian 
caused by loops of virtual new particles (e.g., the last term in ()7.H) ). 

If we find such a deviation from predictions of pure QED (with only electrons and photons), 
we can estimate the scale of "new physics" , i.e. the muon mass M. If there is an unknown 
heavy charged particle, then it interacts with photon and gives a contribution c ~ 
(Fig.EZI)^. 




Figure 57: The muon contribution to the electron magnetic moment 

Operators of the full theory are also expansions in 1/M, in terms of all operators of the 
effective theory with appropriate quantum numbers. In particular, the bare electron and 
the photon fields of the full theory are, up to 1/M^ corrections, 

^o={Q"'i^'o, A,= {CY"A,. (7.2) 
The bare parameters in the Lagrangians of the two theories are related by 

eo={CY'"e',. ao = CX. (7-3) 

^We suppose that it does not interact with electrons directly. 
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The MS renormalized fields and parameters are related by 



(7.4) 



where 

CM 



Z'^{a'{fi),a'{fi)) 



C. aw = #^cs. cM = #^^c:. (7.5) 



The gluon propagators in the two theories are related by 



p2 



+ do 



2\2 



AO 



p2 



+ a. 



°(p2)2 



(7.6) 

This explains why the same decoupling constant C,a describes decoupling for both the 
photon field A and the gauge-fixing parameter a. It is most convenient to do matching at 
— > 0, then the power-suppressed terms in (j7.6j) play no role. The full-theory propagator 
near the mass shell is 

1 



70s 



r - i - n(o) ■ ^^'^"^ 

Only diagrams with muon loops contribute to n(0), all the other diagrams contain no scale. 
In the effective theory 



rylos 



p^ ■' i-n'(o) 

because all diagrams for 11' (0) vanish. Therefore, 

1 



AO 
SA 



At one loop, from ()3.42|1 . 

ZAia) 



1 - 



^A 
^A 



a 
34^ 



n(o)- 



(7.9) 



Z\[a) 



4 a 



3 47r£ 



(because there are two lepton flavours in the full theory). With this accuracy, we may put 
q;'(/x) = a(/x). Re-expressing 11(0) ()6.7|) via a(/i) (and replacing mo by the on-shell muon 
mass M, because the difference is beyond our accuracy), we obtain 



a(/i) 



Z^(a(/i)) 1 - n(o) 



^ _ 4 a(/i) 
3 47r 



M) 



2e 



and flnally 



(7.10) 
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We can always find the /i-dependence from tlie RG equation, tlie initial condition is 
(a{M) = 1 at one loop. 

Let's find Ca{M) with two-loop accuracy. We express Ca ()7.9|) via the renormalized 
quantities: a{M) and the on-shell muon mass M. Technically, this is nearly the same 
calculation as in Sect. 16. 71 But now we express Cg in the one- loop term via a(M), and 
Za{a) contains two lepton flavours: 

(the extra term 16/9 as compared to ()6.62|) comes from the doubled /3o). We can neglect 
the difference between a{M) and a'{M) with our accuracy (this difference is 0{a^), see 
Sect. ESI). From and (lOD we have 



ZA{a) ^ 3 47r£ ^ V 9 ^"^^J \Ane 



All 1/e poles cancel in the renormalized decoupling constant (a{M) ()7.5j) . as they should, 
and we obtain ^ 

a(M) = 1 - 15 f +••• (7-12) 



The RG equation 

dlogfi 



+ 7A(a(/x))-7A(«'(/"))=0 (7.13) 



can be used to find Ca{^^) for /i ^ M. For example, for /i = M(M), the MS muon 
mass (|6.3ip normalized at = M, the difference 

M(M) - M = -4^M , (7.14) 

4:71 



SO that Ca(M(M)) is still 1 + ^(a^): 



13 fa{M) 



2 



a(M(M)) = 1 - - (^^J +••• (7.15) 

This result can be easily obtained directly: if we use Mq = Zm{<y{M))M{M) to renormalize 
the mass in the one-loop term of instead of Mq = Z^M, we get the formula similar 
to (mH), but with a(M(M)) and 13/3 instead of 15. 

7.2 Electron field 

The electron propagators in the full theory and in the low-energy theory are related by 
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It is most convenient to do matching at p 0, where power corrections play no role. The 
full-theory propagator near the mass shell is 

Only diagrams with muon loops contribute to Sv(0), all the other diagrams contain no 
scale; such diagrams first appear at two loops (Fig. EHl)- In the effective theory 

Z'°^ 1 

S-iP)-f, ^r-T3sr(5)-i. ("8) 

because all diagrams for Sv'(O) vanish. Therefore, 

= = (7 IQ) 

z;- 1 - Sv^(o) ■ ^ ■ ^ 




Figure 58: Decoupling of electron field 
At two loops (Fig.EHl), 



{2n 



+ 



(k + p) 



r«eo7 



k^ 



iik g^^ - k^ky)Ii{k 



(7.20) 



where i{k'^g^y—k^kv)T^{k'^) is the muon-loop contribution to the photon self-energy (Fig.l48|l. 
It is transverse; therefore, longitudinal parts of the photon propagators (~ C^a/c/^) do not 
contribute, and the result is gauge invariant. We only need the linear term in p in both 
sides: 



n(fc^) 

{k^Y 



-ten 



d'^k 

(27r)" 



[^^{e^-2p-k^)-f>^ 



+ 2p- k 



Averaging over k directions hj p ■ k jk ^ (k'^/d)^, we obtain 



y ^Ad-l){d-4) f d^k n(P^ 



d 



{2TrY (A;2)2 



(7.21) 
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The muon-loop contribution to the photon self-energy (Fig. I48|) is (we set M = 1; the 
power of M will be restored by dimensionality) 



Contracting in and z/, we obtain 

el f d'^k N 



D^ = l-{k + p)\ D2 = l-k\ N = Tt^^{^ + ^+1)Y{^+1] 
Let's calculate the integral 

d'^p n(p2) 



{27tY (-p2)2 ' 

which appears in ()7.21|) . It is a two-loop massive vacuum diagram of Fig. 

. el f d^ki d'^k2 N 
1 = —I 



(7.22) 



d-lj {2tiY {2tiY D1D2DI ' 
D^ = l-kl, D2 = l-kl, D^ = -{k^-k2)\ 
iV = Tr7^(^i + 1)7^(1^2 + 1). 

All scalar products in the numerator can be expressed via the denominators: 
kl = l-D^, kl = l-D2, k,-k2 = ^{2 + Ds-Di-D2). 

Calculating the trace and omitting terms with Di or D2 (which produce vanishing inte- 
grals), we obtain 

N^2[A-{d- 2)^)3] . 
Our integral / becomes (see ()(i.5(ip ) 

^ = ^(Jp[4ni,l,3)-(ci-2)ni,l,2)]. 

Using ()6.59p and restoring the power of M, we finally arrive at 

_d^m^ = AM-'\ 2(. 2(rf-6) 

{27rY {k^Y (47r)'^ ^ \d - 2){d - 5){d - 7) ' ^' ' 

Therefore, we obtain from (j7.19|) 

.0 _ . ^ etM-'^ 2,. '^{d-l){d-A){d-Q) 

+ {At^Y d{d-2){d-hYd-l)- ' 
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The renormalized decoupling coefficient is 

Its /^-dependence can always be found by solving the RG equation. It is sufficient to obtain 
it at one point, at some specific fi ~ M, to have the initial condition. The most convenient 
point is /i = M, because a{M) = a'{M) + 0{a^) (Sect. ESI) and a(M) = a'(M) + 0{a^) 
(Sect. 17. and the differences can be neglected with our accuracy. The renormalization 
constant Z^ is given by ()5.55|) . and Z[p — by a similar formula with primed coefficients. 
Their ratio is 

Z^{a,a) If 1 „ \ / a 



Z'^{a,a) 



1 + i (^7^0 A/5o + ^ A7^o7;o« - ^7^1^) (4^) ' , (7-26) 



where 7^q is the coefficient of a in 7^,0, and 

4 8 
A/?o = -g , A7^o = 3 , A7^i = -4 

are the single-fiavour contributions to Po, 7^0, Itpi- We obtain 

Z^{a,a) _ / a \2 



) 



Re-expressing ()7.24|) via the renormalized a{M), 



C^ = l + £(l--£ + 



5 \ / a \ 



6 / KAne 



we finally obtain 

Z'^{a,a) hfa{M)V 



The RG equation 

^^^f^ + 7^(«(/i), a(/i)) - 7;(«'(/i), a'(/x)) = (7.28) 
a log jjL ^ 

can be used to find Cv'l/^) 7^ contrast to the case of Ca{,IJ^) ()7.13|) . now 7^ — 7^ 

is of order a^, so that changes of /i or order a (such as, e.g., ()7.14|) ) don't change the 
coefficient of in ()7.27|) . 
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7.3 Charge 

The proper vertex eoT with the external propagators attached (two electron propagators 
5* and one photon propagator D) is the Green function of the fields ipo, ipo, Aq (i.e., the 
Fourier transform of the vacuum average of the T-product of these three fields). Therefore, 
the relation between this quantity in the full theory and in the low-energy effective theory 
is 

eoTSSD = Cj (Ca) e'.T'S'S'D' , (7.29) 
or, taking into account S = C^S', D = C'^D', 

eoT' = iO'' iCir'^" e',T'^ (7.30) 

In the full theory, the vertex at p = p' on the mass shell (p^ = 0) is 

= Z^'Y ■ (7.31) 

Only diagrams with muon loops contribute to K^{p,p)^ all the other diagrams contain no 
scale; such diagrams first appear at two loops (see Fig. l57j) . In the effective theory 

r'^ = , (7.32) 

because all diagrams for A'^{p,p) vanish. Therefore, 

Y>^ = (;Ot"^^ C^^ = ^ = ^r% (7.33) 

and we obtain from (j7.30p 

The situation in QED is simpler, due to the Ward identity fSect. lOj) . In the full theory, 
we have ()6.37|1 Z'^Z^ = 1; similarly, in the effective theory, Z'^^Z^^ = 1 (in fact, these 
two renormalization constants are equal to 1 separately). Therefore, 

CM = 1, C={Cr'. (7.35) 
Recalling also (j3.77j) Z^ = Z~^ and Z'^ = Z'^^, we obtain 

C = C^=fc#V = G^. (7.36) 

Finally, from (I7?T^ . 

UM) = 1 + 15 1- 

47r 

This means that the running charges in the two theories are related by 

a{M) = Ca(M)a'(M) , 



Za V Za 



C(M) = l + 15f4^V + --- (7.37) 
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i.e., the running charge in the full QED (with both electrons and muons) at /i = M is 
slightly larger than in the low-energy effective QED (with only electrons). If you prefer to 
do the matching at /i = M(M) instead of /i = M, you should use 



3 V 47r 

(see dnni)). 



7.4 Electron mass 

In the previous Sections, we considered QED with massless electrons and heavy muons 
(with mass M). Now let's take the electron mass into account as a small correction. We 
shall expand everything up to linear terms in m. The electron propagator in the full theory 
is given by (j3.98j) : in the low-energy effective theory, it involves Sy^ instead of T^y^s- These 



two propagators are related by Cw' 



Ci ^/ /^■?.\ 1 , ^> f^2\ ■ (7.39) 



Comparing the overall factors, we recover ()7.19j) . Comparing the denominators, we obtain 
so that 

mo l-T.y{p^)l + T.'s{p') . Q.-i 1 
''"^ m', l-SV(p2)l + E5(p2) ^^'^^ 1 + S5(0)- ^ 

This is because ^'^.(O) = in our approximation: after we single out tjiq in front of S^(0) 
in the electron self-energy S'(p^) (similar to ()3.97|l ). we can set tjiq = in the calculation 
of ^^(O), and it has no scale. Only diagrams with muon loops contribute to ^^(O); such 
diagrams first appear at two loops fFig. l3Hj). The renormalized mass decoupling constant 
is 

m'(/i) Zm{a{^)) 

At two loops (Fig. EHI), we have to take mo into account, as compared to the massless 
case (Hm): 

-'^^P^ = J J2^^'''^ \k + py-mf ''^ ) ' ^^^^^^^^ ^ ■ 

Again, this is gauge-invariant. In order to extract S^, we have to retain mo in the numerator 
of the electron propagator; after that, we can put mo = 0: 



T\2 
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Using the integral ()7.23|) . we obtain 



y(0)- ^o^"'V 2. . 2{d-l){d-Q) 

Therefore, from (17^ and ^HM . 

. , 8(d-l)(d-Q) 

(^0 = 1 + _o Y^e) ^ f7 44) 

^"^ [Any ^U{d-2){d-5){d-7) ^ ' 

It is most easy to calculate Cm(^); Cmilj) can then be found by solving the RG equation. 
To our accuracy, we can take a{M) = a'{M); the ratio of the renormalization constants 
Zm{c() / Z!^{a) is given by the formula similar to ()7.26p . but •y^ is gauge-invariant (7^q = 0), 
and (see (f^.^Sjll 7^0 = 6, A-fmi = -20/3: 

Zm(ce) ^ ( ^ 5 \ / a 



2--e 



ZL(a) V 3 / KAne 



Re-expressing (j7.44p via the renormalized a{M), 



we finally obtain 



The RG equation 



18 / V47re 



C„(M) = 1- — M-^ +■■■ (7.45) 



18 V A-n 

rflogCm(/u) 



+ 7,n(a(/i))-7™(a'(/i))=0 (7.46) 



dlogyU 

can be used to find Cm(/^) for yU 7^ M (the difference 7.^ — 7^ is of order a^). 



7.5 Decoupling in QCD 

In QED, effects of decoupling of muon loops are tiny. Also, pion pairs become important at 
about the same energies as muon pairs, so that QED with electrons and muons is a model 
with a narrow region of applicability. Therefore, everything we discussed in the previous 
Sections is not particularly important, from the practical point of view. 

In QCD, decoupling of heavy flavours is fundamental and omnipresent. It would be 
a huge mistake to use the full 6-flavour QCD at characteristic energies of a few GeV, or 
a few tens of GeV: running of ^^(/i) and other quantities would be grossly inadequate, 
convergence of perturbative series would be awful because of huge logarithms. In most 
cases, anybody working in QCD uses an effective low-energy QCD, where a few heaviest 
flavours have been removed. Therefore, it is important to understand decoupling in QCD. 
And to this end the lessons of QED are very helpful. 
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Suppose we have a heavy flavour with on-shell mass M and rii hght flavours. Then 
running of the full-theory coupling a^J^^^^\ij) is governed by the {rii + l)-flavour /3-function; 
running of the effective-theory coupling ai^'\fi) is governed by the ra^-flavour /3-function; 
their matching is given by 

= , (7.47) 

with 

32^ \ fas{M)^' 



C^{M) = l+\^15CF--CAjTF\^-^^j +•■■ (7.48) 

Here the Cp term can be obtained from the QED result ()7.37p by inserting the obvious 
colour factors; the Ca term is more difficult to obtain. The RG equation 

"^^"f^"^^^ + 2/5("'+i)(«i"'"''H/^)) - 2p^^^\ai^^\fi)) = (7.49) 
ct log /i 

can be used to find CaifJ') for /i ^ M. The difference 

and therefore for the MS mass M(M) (jfi.fiC'il) we have the well-known formula [2H! 

C„(M(M)) = 1 + (He, - T, {^)\ ■ ■ ■ (7,50) 

The QCD running coupling ^^(yu) not only runs when fi varies; it also jumps when 
crossing heavy-flavour thresholds. The behaviour of as(/i) near rub is shown in Fig. EHl 
(this figure has been obtained using the Mathematica package RunDec [21], which takes 
into account 4- loop /5-functions and 3-loop decoupling). At fi > rrih, the correct theory is 
the full 5-flavour QCD (ai^-'(/i), the solid line); at /i < mi,, the correct theory is the effective 
low-energy 4-flavour QCD {ai^\fi), the sohd line); the jump at /i = (I7.48p is shown. Of 
course, both curves can be continued across nib (dashed lines), and it is inessential at which 
particular /i ~ rub we switch from one theory to the other one. However, the on-shell mass 
nib (or any other mass which differs from it by 0{as), such as, e.g., the MS mass nib{nib)) 
is most convenient, because the jump is small, 0{a^). For, say, fj, = 2nib or fi = nib/2 it 
would be 0{al). 

Light-quark masses nii{fi) also rum with /i, and also jump when crossing a heavy-quark 
threshold. The QCD result 



m("'+^)(M) =m("')(M) 



89 fas{M) 

— Lyp-^F 



18 V 47r 



(7.51) 



can be obtained from the QED one ()7.45|) by inserting the obvious colour factors. 
All QCD decoupling relations are currently known with three-loop accuracy 0. 
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0.225 n 




0.21 -I > , 

nib — O.SGeV m^, + 0.5GeV 



Figure 59: as {^) and as 

8 Conclusion: Effective field theories 

We already discussed effective field theories in Sect. When using QCD, one rarely works 
works with all 6 flavours; more often, one works in an effective low-energy QCD with fewer 
flavours. 

In fact, all our theories (except The Theory of Everything, if such a thing exists) 
are effective low-energy theories. We don't know physics at arbitrarily small distances 
(maybe, even our concept of space-time becomes inapplicable at very small distances). We 
want to describe phenomena at distances larger than some boundary scale; our ignorance 
about very small distances is parametrized by local interactions of low-mass particles. 
These observable particles and their interactions are thus described by an effective field 
theory with all possible local operators in its Lagrangian. Coefficients of higher- dimensional 
operators have negative dimensionalities, and are proportional to negative powers of the 
energy scale of a new physics (the scale at which the effective low-energy theory breaks 
down). At energies much lower than this scale, these higher-dimensional terms in the 
Lagrangian are unimportant. We may retain dimension-4 terms (renormalizable) and, 
maybe, one or two power corrections. 

The first historical example of an effective low-energy theory is the Heisenberg-Euler 
effective theory in QED. It is still the best example illustrating typical features of such 
theories. 

In order to understand it better, let's imagine a country, Photonia, in which physicists 
have high-intensity sources and excellent detectors of low-energy photons, but they don't 
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have electrons and don't know that such a particle exists^. At first their experiments 
(Fig. lUUk ) show that photons do not interact with each other. They construct a theory, 
Quantum Photodynamics, with the Lagrangian 

L = -iF^,F'^^ (8.1) 

But later, after they increased the luminosity (and energy) of their "photon colliders" and 
the sensitivity of their detectors, they discover that photons do scatter, though with a very 
small cross-section fFig.l60b). They need to add some interaction terms to this Lagrangian. 
Lowest-dimensional operators having all the necessary symmetries contain four factors F^y. 
There are two such terms: 

L = -\f^.F^^' + ci {F^yF^'f + C2F^yF''"F^fsF'"^ . (8.2) 

They can extract the two parameters ci,2 from two experimental results, and predict results 
of infinitely many measurements. So, this effective field theory has predictive power. 




a b 



Figure 60: Scattering of low-energy photons 

We know the underlying more fundamental theory for this effective low-energy the- 
ory, namely QED, and so we can help theoreticians from Photonia. The amplitude of 
photon-photon scattering in QED at low energies must be reproduced by the effective 
Lagrangian ()8.2|) . At one loop, it is given by the diagram in Fig. |^ Expanding it in the 
photon momenta, we can easily reduce it to the massive vacuum integrals ()2.1|) . Due to 
the gauge invariance, the leading term is linear in each of the four photon momenta. Then 
we equate this full-theory amplitude with the effective-theory one following from (|8.2p , and 
find the coefficients ci,2 (this procedure is known as matching). The result is 

L = -\f^.F^^ + [-5 {F,yFn^ + lAF^yF'^^F^pF^^] . (8.3) 

It is not (very) difficult to calculate two-loop corrections to this QED amphtude using the 
results of Sect. 16.61 and thus to obtain terms in these coefficients. 

^We indignantly refuse to discuss the question "Of what the experimentahsts and their apparata are 
made?" as irrelevant. 
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Figure 61: Photon-photon scattering in QED at one loop 



There are many apphcations of the Lagrangian ()8.3p . For example, the energy density 
of the photon gas at temperature T is ~ by dimensionality (Stefan-Boltzmann law). 
What is the radiative correction to this law? Calculating the vacuum diagram in Fig. 1^21 
at temperature T, one can obtain [HD] a correction ~ a'^T^/m*. Of course, this result is 
only valid at T ^ m. 



The interaction terms in the Lagrangian (j8.3|) contain the "new physics" energy scale, 
namely the electron mass m, in the denominator. If we want to reproduce more terms in 
the expansion of QED amplitudes in the ratio u/m {uj is the characteristic energy), we 
can include operators of higher dimensions in the effective Lagrangian; their coefficients 
contain higher powers of m in the denominator. Such operators contain more F^,^ and/or 
its derivatives. Heisenberg and Euler derived the effective Lagrangian for constant field 
containing all powers of F^j^; it is not sufficient for finding coefficients of operators with 
derivatives of F^^,. The expansion in u/m breaks down when ~ m. At such energies 
the effective low-energy becomes useless, and a more fundamental theory, QED, should be 
used; in particular, real electron-positron pairs can be produced. 

QED, the theory of electrons and photons, is an effective low-energy theory too. The 
leading non-renormalizable term in its Lagrangian has dimensionality 5. It is the Pauli 
interaction, the last term in (|7.1|) . It contains the scale of new physics, in this case the 
muon mass M, in the denominator. In reality, hadronic loops produce a comparable contri- 
bution to the electron magnetic moment. This hadronic contribution cannot be calculated 
theoretically, because no one knows how to calculate low-energy hadronic processes from 
the first principles of QCD^. Therefore, the coefficient of the Pauli term in the QED La- 
grangian is a phenomenological parameter, to be extracted from experiment — from the 
measurement of the electron magnetic moment. 

We were lucky that the scale of new physics in QED is far away from the electron 

^It can be calculated, with a good accuracy, from experimental data about e^e~ annihilation into 
hadrons. 




Figure 62: Radiative correction to the Stefan-Boltzmann law 
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mass m. Contributions of heavy-particle loops are also strongly suppressed by powers of a. 
Therefore, the prediction for the electron magnetic moment from the pure QED Lagrangian 
(without non-renormalizable corrections) is in good agreement with experiment. After this 
spectacular success of the simplest Dirac equation (without the Pauli term) for electrons, 
physicists expected that the same holds for the proton, and its magnetic moment is e/ (2mp). 
No luck here. This shows that the picture of the proton as a point-hke structureless particle 
is a very poor approximation at the energy scale rrip. 

Another classical example of low-energy effective theories is the four-fermion theory 
of weak interactions. It was first proposed by Fermi as a fundamental theory of weak 
interactions; now it is widely used (in a shghtly modified form) as an effective theory of 
weak interactions at low energies. 

Let's consider, for example, 6-quark decays. In the Standard Model, they are described 
by diagrams with W exchange. However, the characteristic energy E ~ mf, ^ myy, and 
we can neglect the W momentum as compared to mw, and replace its propagator by the 
constant In other words, the W propagation distance ~ 1/mw is small compared 

to characteristic distances of the process ~ We can replace the W exchange by a 

local four-fcrmion interaction; it is described by dimension-6 four-fermion operators, their 
coefficients contain There are 10 such operators. We can calculate QCD corrections 

of any order to diagrams containing such four-fermion vertex without any problems; UV 
divergences are eliminated by the matrix of renormalization constants of these operators. 
Their coefficients in the effective low-energy Lagrangian obey a renormalization group 
equation. The initial condition at /i — mw is obtained by matching the full-theory results 
with the effective-theory ones. 

The "new physics" scale where this effective theory breaks down is mw As is always 
the case in effective theories, we have to fix the order of the expansion in E'^/m'^ we are 
interested in. In this case, we work at the first order in and there can be only one 

four-fermion vertex in any diagram. What if we want to include the next correction in 
E'^/m^rl In this case we have to include all the relevant dimension-8 operators into the 
Lagrangian; their coefficients contain Their number is large but finite. Elimination 

of UV divergences of diagrams with two four-fermion vertices (i.e., renormalization of 
products of dimension-6 operators) requires not only renormalization of these operators 
separately, but also local overall counterterms (i.e., dimcnsion-8 operators). 

We see that the fact that the theory with four-fermion interaction is not renormalizable 
(as all effective field theories) does not prevent us from calculating finite renormalized 
results at any fixed order in the expansion in E"^ /m^. The number of coefficients in the 
Lagrangian is finite at each order of this expansion, and the theory retains its predictive 
power. In practice, this number quickly becomes large, and one can only deal with few 
power corrections in an effective theory. 

In the past, only renormalizable theories were considered well-defined: they contain a 
finite number of parameters, which can be extracted from a finite number of experimental 
results and used to predict an infinite number of other potential measurements. Non- 
renormalizable theories were rejected because their renormalization at all orders in non- 
renormalizable interactions involve infinitely many parameters, so that such a theory has 
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no predictive power. This principle is absolutely correct, if we are impudent enough to 
pretend that our theory describes the Nature up to arbitrarily high energies (or arbitrarily 
small distances). 

Our current point of view is more modest. We accept the fact that our theories only 
describe the Nature at sufficiently low energies (or sufficiently large distances). They are 
effective low-energy theories. Such theories contain all operators (allowed by the relevant 
symmetries) in their Lagrangians. They arc necessarily non-renormalizable. This does not 
prevent us from obtaining definite predictions at any fixed order in the expansion in E /M, 
where E is the characteristic energy and M is the scale of new physics. Only if we are lucky 
and M is many orders of magnitude larger than the energies we are interested in, we can 
neglect higher-dimensional operators in the Lagrangian and work with a renormalizable 
theory. 

Practically all physicists believe that the Standard Model is also a low-energy effec- 
tive theory. But we don't know what is a more fundamental theory whose low-energy 
approximation is the Standard Model. Maybe, it is some supersymmetric theory (with 
broken supersymmetry) ; maybe, it is not a field theory, but a theory of extended objects 
(superstrings, branes); maybe, this more fundamental theory lives in a higher-dimensional 
space, with some dimensions compactified; or maybe it is something we cannot imagine at 
present. The future will tell. 

I am grateful to K.G. Chetyrkin and V.A. Smirnov for numerous discussions of the 
topics discussed in these lectures; to the organizers of the DIAS TH-2005 school in Dubna 
for inviting me to give this course; to S.V. Mikhailov and A.A. Vladimirov for useful 
remarks on the contents of the lectures. 
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A Colour factors 

Here we shall discuss how to calculate colour diagrams for SU{Nc) colour group (in the 
Nature, quarks have Nc = 3 colours). For more details, see 

Elements of SU{Nc) are complex A^^^ x A^^ matrices U which are unitary {U^U = 1) 
and have detU = 1. Complex A^^^- component column vectors g* transforming as 

q^Uq or q'-^U'jq^ (A.l) 

form the space in which the fundamental representation operates. Complex conjugate row 
vectors q^ = (g*)* transform as 

q+^q+U+ or qt^qfiU+Yi where {U+Y i = {U' j)* ; (A.2) 

this is the conjugated fundamental representation. The scalar product is invariant: q^q' 
q~^U^Uq' = q~^q'. In other words, 

5] ^ 5ff/\(t/+)', = f/\(f/+)^- = Si (A.3) 

is an invariant tensor (it is the colour structure of a meson). The unit antisymmetric 
tensors e*i---*^c and ei-^^,,,i^^ are also invariant: 

gii-iivc ^ ^ii-iive u''.^ . . . U'^'-j^^ =detU- e''-'''- = e'^-'''- (A.4) 

(they are the colour structures of baryons and antibaryons). 
Infinitesimal transformations are given by 

f/ = l+ia''r, (A.5) 

where a°' are infinitesimal real parameters, and t"" are called generators (of the fundamental 
representation). They have the following properties: 

f/+f/ = 1 + ia'^ (t'* - (r)+) = 1 ^ (t'^)+ = r 

(A.6) 

det f/ = 1 + ia'' Tr t"^ = 1 ^ Tr = , 

and are normalized by 

Tr t'^t'' = Tp^'^^ (A. 7) 

usually, Tp = 1/2 is used, but we shall not specialize it. The space of unitary matrices is 
A'^-dimensional, and that of traceless unitary matrices — (A^^ — l)-dimensional. Therefore, 
there are A"^ — 1 generators which form a basis of this space. Their commutators are i 
times unitary traceless matrices, therefore, 

[e,t^]=iP^H' , (A.8) 

where /"^^^ are real constants. 
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The quantities 
transform as 

A" q+U+fUq' = W^A^ ; (A.9) 
this is the adjoint representation. It is defined by 

U+t''U = U''H\ (A.IO) 

and hence 

^ab = ^ Tr U+eUt^ . (A. 11) 
Tf 

The components are some fixed numbers; in other words, they form an invariant 

tensor (see ()A.10|) ): 

^ V'^Wkit^fiiU+y J = {fYj . (A.12) 
For an infinitesimal transformation, 

A"" g+(l - ia^t^)t%l + ia^t'')q' = q+{t'' + ia'ir'h'')q' , 

so that 

jjab ^ gab ^ -^c^^cyb ^ (A.13) 

where the generators in the adjoint representation are 

As for any representation, these generators satisfy the commutation relation 
it follows from the Jacobi identity 

It is very convenient to do colour calculations in graphical form 0]. Quark lines mean 
(5*, gluon lines mean 5°^, and quark-gluon vertices mean (t")*j. There is no need to invent 
names for indices; it is much easier to see which indices are contracted — they are connected 
by a line. Here are the properties of the generators t"^ which we already know: 

Trl = iV, or =N,, 

TrT = or 'H^^^ = ° ' {P-AQ) 

TitH^ = Tf5''^ or -w^^^^^^ = Tp -wwv, . 
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There is a simple and systematic method for calculation of colour factors — Cvitanovic 
algorithm [l]. Now we are going to derive its main identity. The tensor {t°'yj{t°'Yi is 
invariant, because is invariant. It can be expressed via (5*-, the only independent 

invariant tensor with fundamental-representation indices (it is clear that and eii...jjvc 

cannot appear in this expression, except the case Nc = 2; in this case e'^'^eji can appear, 
but it is expressible via 5*). The general form of this expression is 

{n,{n'i = a[6l6^-b6}6^] , (A.17) 

or graphically 

> ; > — " — 



< > < 



- b 



where a and b are some unknown coefficients. If we multiply ()A.17jl by S: 

iniin'i = = a[6f- bNjf] , 
i. e., close the upper line in flA.18|l . 



< < < 



6 



we obtain 



If we multiply (IXTfll by (t^i, 

i. e., close the upper line in (lA.lSjl onto a gluon, 

1 



1 




we obtain 



The final result is 



(A.18) 



(A.19) 
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or graphically 



1 



(A.20) 



This identity allows one to eliminate a gluon exchange in a colour diagram: such an 
exchange is replaced by the exchange of a quark-antiquark pair, from which its colour- 
singlet part is subtracted. 
It can be also rewritten as 



or 



1 




1 

7> 



(A.21) 



(A.22) 



This shows that the product of the fundamental representation and its conjugate is re- 
ducible: it reduces to the sum of two irreducible ones, the adjoint representation and the 
trivial one. In other words, the state of a quark-antiquark pair with some fixed colours is 
a superposition of the colour- singlet and the colour-adjoint states. 

The Cvitanovic algorithm consists of elimination gluon exchanges ()A.20|) and using 
simple rules fjA.16|) . Let's consider a simple application: counting gluon colours. Their 
number is 




as we already know. 

Now we consider a very important example: 




The result is 



or rr = CF: 



(A.23) 
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where the Casimir operator in the fundamental representation is 

Cf = Tf (^N, - . (A.24) 

Colour diagrams can contain one more kind of elements: 3-gluon vertices if"'^'^. The 
definition ()A.8jl when written graphically is 




(A.25) 



Let's close the quark hne onto a gluon: 



Therefore, 




(A.26) 



This is the final rule of the Cvitanovic algorithm: elimination of 3-gluon vertices. 
The commutation relation ()A.15|) can be rewritten graphically, similarly to ()A.25|) : 




(A.27) 



Sometimes it is easier to use this relation than to follow the Cvitanovic algorithm faithfully. 
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Now we consider another very important example: 




-'-F 



2 





2 





— 2TpNc '^y^/^y^ . 



The result is 




or iP'^Hf'^'' = Ca5''\ 



where the Casimir operator in the adjoint representation is 
Another example: 



> «.> i> »> ^Tw 




1 



> <> 



Tf I 
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One more example: 




(A.32) 



> > > — 

or 

if^'HY = . (A.33) 
This can be derived in a shorter way. Using antisymmetry of /°*'^, ()A.25|) and ()A.28|) . 




(A.34) 

We shall also need the result 




(A.35) 



It can be, of course, obtained by using the Cvitanovic algorithm. A shorter derivation is 
similar to ()A.34|) . with the lower quark line replaced by the gluon one; the commutation 
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relation ()A.27|) in the adjoint representation is used instead of ()A.25|) in the fundamental 
one. 

Sometimes the colour factor of a diagram vanishes, and we don't need to undertake its 
(maybe, difficult) calculation. For example. 




(A.36) 



This can be obtained by using the Cvitanovic algorithm. A shorter derivation is following. 
There is only one colour structure of the quark-gluon vertex: 




Therefore, let's close this diagram: 



; ^^^^^^^^^^^ 






Colour factors of not too complicated diagrams (e.g., propagator and vertex diagrams 
which have to be calculated for obtaining QCD fields anomalous dimensions and /3-function) 
with L loops (including Lq quark loops) have the form 



(A.37) 



This form is valid not only for SU{Nc), but for any gauge group. At Lg > 4, new Casimir 
operators appear, which are not reducible to Cp and Ca- Suppose we have calculated the 
colour factor for SU {Nc) using the Cvitanovic algorithm, and we want to rewrite it in the 
form ()A.37|) valid for any group. Our colour factor is (TpUj)^'' times a sum of terms of the 
form Tp'N^ with various values of n. Using ()A.24j) and ()A.29j) . we can rewrite each such 
term as ^ ^ ^ 

^) ~ '^^) ^^^^^ x + y = Lg 

Therefore, we may substitute 



X 



y 



n . 



Ca 
2 



1 - 



2Cf 



2Cf 



(A.38) 



Of course, if there are no errors, the result will contain no fractional powers, just a poly- 
nomial in and Ca-^^ 



am grateful to K.G. Chetyrkin for explaining this trick to me. 
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